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I.  INTRODUCTION 


With  the  advent  of  the  missile  era  has  come  a  revolution 
in  naval  warfare  methods.  The  problem  of  contending  with 
"smart"  weapons  has  given  a  whole  new  meaning  to  both 
offense  and  defense  alike.  These  missiles  may  employ  a 
variety  of  acquisition  hardware.  Some  will  be  passive,  some 
active,  and  some  both. 

Particular  problems  arise  when  missiles  use  passive 
acquisition  methods.  Passive  sensors  radiate  nothing  and 
are  therefore  invisible  to  conventional  electronic  support 
measures  (ESM) .  Furthermore,  their  small  size  allows  them 
to  reflect  only  meager  amounts  of  radar  energy,  which 
further  reduces  reaction  and  countermeasure  employment  time 
of  the  ship.  The  countermeasures  may  be  in  the  form  of 
onboard  or  off board  weapons  systems  or  deception  devices. 

The  cost  versus  effectiveness  of  onboard  close-in  weapon 
systems  has  been  a  driving  factor  in  the  Navy's  desire  to 
analyze  the  employment  of  off board  expendable  seduction  and 
distraction  devices  as  a  relatively  inexpensive  alternative. 

Here  we  will  restrict  our  studies  to  the  passive 
infrared  seeking  missile  problem,  although  concurrent 
analysis  is  being  conducted  to  contend  with  the  active  RF 


missile  seekers 


To  reduce  radically  the  cost  o-f  development  of  off  board  • 
device  characteri sties  and  deployment  tactics,  a  computer 
simulation  has  been  constructed  to  assist  in  the  analysis 
process.  This  modelling  program  is  appropriately  named  the 
“Stochastic  Infrared  Engagement  Model  -  SIREM" .  It 
originated  at  Naval  Research  Laboratory  in  Washington,  D.C. , 
and  was  transported  to  the  Naval  Postgraduate  School  in 
Monterey,  CA  to  assist  in  the  analysis  of  the  use  of 
infrared  decoys  as  a  counter  to  the  passive  seeker  problem. 

It  is  designed  to  imitate  the  typical  cruise-missile  versus 
ship  engagement  scenario  and  to  establish  measures  of 
effectiveness  (MOEs)  for  both  1)  decoy  characteristics,  and 
2)  deployment  tactics.  Therein  lies  the  basis  for 
continuing  verification  and  upgrades  to  the  model  so  that  it 
remains  capable  of  accurate  evaluation  of  the  hardware  at 
hand. 

Most  real-world  systems  are  very  complex.  In  order  to 
model  any  system  accurately,  all  factors  which  affect  it 
must  be  represented  in  as  accurate  a  fashion  as  possible. 
Weaknesses  in  representati on  should  be  identified  and 
documented  so  that  overconfidences  in  the  output  do  not 
occur.  Once  documented,  the  weak  areas  of  the  model  may  be 
studied  and  reworked  to  increase  the  accuracy  of  the  model. 

This  paper  deals  with  the  evolution  of  SIREM.  By 
simulating  each  component  of  the  engagement  scenario  in  a 
realistic  manner,  the  optimum  decoy  deployment  schemes  can 


and 


be  obtained,  combined  with  their  RF  counterparts, 
algorithmically  automated  -for  deployment  via  computer  in 
order  to  reduce  chances  o-f  operator  control  errors.  It  can 
easily  be  seen  here  that  it  is  not  financially  or 
logistical ly  feasible  to  evaluate  all  scenarios  in  field 
tests  using  real  ships  and  missiles,  since  the  cost  Mould  be 
prohibitive. 

Extensive  conversion  and  review  of  the  SIREM  model  was 
conducted  to  assist  in  its  enhancement.  Two  specific  areas 
of  weakness  were  identified  which  require  some  attention. 

The  first  area  of  concern  was  a  lack  of  atmospheric 
attenuation  on  the  emitted  IR  signatures  of  the  ship  and/or 
decoys.  The  second  regards  the  study  of  seeker  head 
di scri mi nation  techniques  and  how  to  implement  them  based  on 
data  available.  The  emphasis  in  this  paper  revolves  around 
finding  a  simple  technique  for  calculating  atmospheric 
transmittance  throughout  the  path  of  an  incoming  infrared 
seeking  missile  which  will  accurately  affect  its  probability 
of  acquisition  in  the  simulation  environment.  Two  methods 
currently  used  in  obtaining  transmittance  in  a  stochastic 
model  such  as  SIREM  involve  either  look-up  tables  which  are 
either  incomplete  or  inflexible,  or  incorporation  of  the 
LOUTRAN  program  as  a  subroutine.  The  second  method  is 
preferred  since  it  is  very  flexible  to  observed  atmospheric 
conditions;  however,  it  requires  an  extension  of  9000  lines 
of  FORTRAN  code  along  with  the  added  input-output  coding. 
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which  serves  only  to  calculate  one  value  of  atmospheric 
transmittance  for  each  incremental  range  throughout  the 
missile  flight  path.  It  was  felt  by  the  author  that  this 
was  excessively  wasteful  of  valuable  computer  time  and 
memory  requirements  considering  the  model's  current  size  of 
approximately  15,000  lines  of  code. 

The  intent  was  to  devise  a  simpler  scheme  for  finding  as 
accurately  as  possible  what  the  atmospheric  transmittance 
would  be  for  a  given  set  of  atmospheric  conditions  that 
would  be  easily  incorporated  into  the  incremental  flight 
path  of  a  missile  and  more  accurately  affect  the  seeker  head 
acquisition  probability. 

To  more  fully  understand  the  principle  of  atmospheric 
transmi ttance  and  its  effects,  one  chapter  is  dedicated  to 
explaining  the  basic  principles  of  atmospheric 
transmi ttance.  Once  a  good  understanding  of  the  principles 
of  atmospheric  transmittance  has  been  gained  by  the  reader, 
an  introduction  to  the  L OUTRAN  program  is  performed, 
presenting  its  capabilities  for  prediction  of  transmittance. 
This  will  be  followed  by  a  simplistic  approach  for  modeling 
atmospheric  transmittance  with  very  short  computer  coding 
schemes. 

The  closing  chapters  will  be  dedicated  to  more  detailed 
explanation  of  the  SIREN  modelling  program,  followed  by  an 
introduction  to  possible  di scrimination  techniques  Since 
seeker  di scrimination  adds  to  the  problems  surrouding  the 


use  of  off board  countermeasures,  the  subject  of  enhancing 
current  decoy  designs  by  taking  advantage  of  certain  natural 
phenomenon  is  addressed  in  the  concluding  chapter. 

Future  enhancements  to  SIREM  to  incorporate  both 
transmittance  and  discrimination  techniques  are  the  object 
of  this  work  and  concurrent  work  by  Naval  Research 
personnel.  These  enhancements  are  intended  to  more 
accurately  reflect  the  acquisition  probability  of  the 
seekers  being  modelled  so  that  refinements  to  the  design  and 
deployment  of  offboard  countermeasures  are  possible. 


1 1 .  ATMOSPHERIC  ENV I RQNMENT 

This  chapter  is  dedicated  to  reviewing  the  principles  of 
atmospheric  propagation  of  infrared  radiation,  the  causes  of 
its  variability,  and  tools  available  to  predict  how  much  is 
transmitted  at  various  ranges  and  atmospheric  conditions. 

The  atmosphere  is  comprised  of  numerous  gasses  - 
nitrogen,  oxygen,  water  vapor,  carbon  dioxide,  methane, 
carbon  monoxide,  and  ozone  to  name  the  major  components. 
Additionally  there  are  numerous  suspended  particles 
distributed  in  various  densities  throughout  the  earth's 
atmosphere.  The  major  constituents  of  the  atmosphere  are 
nitrogen,  comprising  78%,  and  oxygen,  which  occupies  20%  of 
the  total  volume.  Water  vapor  makes  up  1%,  while  carbon 
dioxide  constitutes  only  0.04%  of  the  atmosphere.  The 
transmission  of  infrared  radiation  through  the  atmosphere 
depends  on  the  meteorological  conditions,  and  hence  varies 
with  weather  conditions  and  altitude.  Only  an  approximate 
result  can  be  gained  through  the  theoretical  study  of  the 
principles  of  atmospheric  attenuation  since  no  definitive 
conclusions  can  be  made.  The  base  of  knowledge  regarding 
atmospheric  attenuation  and  the  associated  formulae  comes 
principally  from  empirical  means  or  approximations. 


The  primary  -factor  in  absorption  in  the  in-frared 
spectrum  is  water  vapor.  The  wide  variations  in 
concentration  and  distribution  o-f  water  vapor  in  the 
atmosphere  are  the  primary  reasons  -for  the  variance  in 
absorption  levels  -from  one  meteorological  profile  to  the 
next.  Although  carbon  dioxide  comprises  only  4/100  percent 
volume  o-f  the  atmosphere,  it  is  second  in  importance  as  an 
attenuator  of  infrared  radiation,  as  well  as  being  more 
uniformly  distributed  than  water  vapor. 

There  are  three  major  factors  which  contribute  to  the 
attenuation  of  infrared  radiation:  Extinction,  or  the  loss 
of  energy  resulting  from  the  interaction  of  the  beam  with 
the  various  absorbing  or  scattering  constituents  contained 
in  the  atmosphere;  Refraction,  or  the  bending  of  light  rays 
due  to  the  refractive  index  gradients  of  the  atmosphere,  and 
Scintillation,  or  the  distortion  of  the  optical  beam  due  to 
small  scale  turbulence.  Refraction  will  not  be  discussed 
heavily  here  since  it  will  be  apparent  later  that  it  may 
serve  only  to  extend  the  optical  horizon  of  infrared  seeking 
missiles,  but  have  little  affect  on  the  simulation  discussed 
herein.  Extinction  represents  the  combined  effects  of  the 
scattering  of  radiation  into  or  out  of  the  beam  due  to 
aerosols  and  air  molecules  present  in  the  beam  path,  and 
absorption  caused  by  the  direct  transfer  of  energy  from 
radiation  generated  by  the  source  to  the  vibrational  or 
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rotational  kinetic  energy  o-f  the  molecules  and  aerosols. 
This  is  illustrated  in  the  -following  diagram. 


Transmitted 
energy, 


4  Energy  scattered  out  o-f  beam 
yC- Energy  scattered  into  beam 

|  IR  |  tt>0  =  »>|  Receiving  I 
I  Source  I  =  =  =  =  =  =  =  =  «>!  IR  I 
I  Radiator  I  =  »>9  0=  «>|  Detector  I 


Absorbed  energy 


Beam 


^  l 


Emitted  energy 


Figure  2-1.  Causes  o-f  extinction 


The  attenuation  or  extinction  of  electromagnetic 
radiation  is  described  by  the  Lamb ert -Beer -Bouguer  Law 


1  =  1.  €  *" 


(Eq.  2.1) 


where  I  is  the  attenuated  radiation,  I.  is  the  source 
radiation,  p  is  the  attenuation  coefficient,  and  R  is  the 
path  length.  The  extinction  is  caused  by  absorption  and 
scattering  by  molecules  and  aerosols.  They  contribute 
linearly  to  the  the  total  extinction  as  follows: 


P  —  PnA  +  Pm  +  P»n  +  Pm  (Eq.  2.2) 

where  Pm-  molecular  absorption  by  water  vapor,  ozone,  etc 
#m=  molecular  (Rayleigh)  scattering 
Pm=  aerosol  absorption  (dry  particles  neglected) 

£••=  aerosol  (die)  scattering 


The  importance  of  each  of  these  effects  depends  on  the 
wavelength  of  the  electrooptic  energy.  The  wavelength 


dependencies  exist  because  o-f  the  relative  sizes  of  the 
responsible  constituents  (Molecule  or  aerosol)  and  the 
energy  states  of  the  Molecules  (Molecular  absorption  bands). 
The  importance  of  the  principal  component  of  attenuation  for 
the  primary  wave  bands  is  given  here. 


TABLE  2.1.  PRIMARY  ATTENUATORS  IN  IR  SPECTRUM 


Wavelength  Region 
(  In  micrometers  ) 

Visible  .4  -  .7  micron 

Near  IR  .7  -  2  micron 

Mid  IR  3-5  micron 

Far  IR  8-14  micron 


Attenuation  Coefficients 
In  order  of  importance 

Pm  i  Pm 

Pm  iPm  for  some  wavelengths) 
Pm  ( HjO )  f  Pm 
Pm  (HjO)  ,  Pm 


The  aerosol  relationship  used  to  describe  extinction  is 
principally  governed  by  the  size  (aerosol  radius)  of  the 
particle  and  the  wavelength  of  the  incident  energy  as: 


1)  Rayleigh  Scattering  (where  r<<X  -  molecule  effects). 

2)  Mie  Scattering  (where  r  »  X  approximately). 

3)  Geometric  Scattering  (where  r>>\  -  large  particles). 

An  illustration  of  the  scattering  description  is  as  shown 
here. 

I(*)sScattered  light 
Scattering  Angle 

I  1  Incident  Energy _ _ <)“0* 


Figure  2-2.  Scattering  Geometry 


the  principal  factor  in  the  sky's  blue  appearance.  Q  is 
higher  at  the  longer  wavelengths.  For  geometric  scattering 
Q  *  2.0.  Only  particles  with  radii  between  about  0.1  and 
1.0  micrometers  are  effective  in  scattering  light  in  the 
visible  range.  These  particles  are  too  small  to  rain  out 
due  to  their  low  colle'cion  efficiency.  This  is  why  rain 
alone  does  not  reduce  haze  and  increase  visibility.  A 
charge  in  the  vertical  temperature  profile  (such  as  the 
breaking  up  of  an  inversion  layer  by  a  frontal  passage)  does 
increase  horizontal  visibility  through  upward  diffusion  of 
the  aerosols. 

Because  the  scattering  efficiency  Q  is  a  function  of 
aerosol  size  and  wavelength,  scattering  itself  is  a  function 
of  wavelength. 

»  X*  (Eq.  2.5) 

Here  a  is  the  Angstrom  coefficient,  and  varies  generally 
from  0.5  to  2.0.  This  value  may  be  as  high  as  4.0  for 
Rayleigh  scattering.  The  scattering  coefficients  for 
various  aerosol  types  are  shown  below. 
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With  this  diagram  we  see  that  the  wavelength  dependence  of 
scattering  explains  why  the  sky  is  blue,  since  the  shorter 
wavelengths  are  scattered  so  much  more.  Likewise  heavy  haze 
can  cause  the  red  sunset  effect  since  the  longer  wavelengths 
are  not  scattered  out  of  the  path.  Note  also  how  the  cloud, 
fog,  and  haze  coefficients  carry  over  into  the  infrared 
region.  Me  can  expect  these  to  be  major  factors  in 
transmittance  reduction  there. 

Scattering  is  also  indirectly  a  function  of  relative 
humidity  (RH)  because  as  RH  increases,  the  water  molecules 
collect  on  a  dry  particle  until  enough  of  them  are  on  the 
particle  to  cause  it  to  go  into  a  solution  and  form  a 
droplet.  This  is  called  deliquescence.  Since  a  solution 
droplet  can  grow  in  size  much  faster  than  a  particle,  the 
scattering  increases  rapidly  with  RH  due  to  an  aerosol 
rapidly  changing  its  size  distribution  through  growth  of 
droplets.  Typical  scattering  coefficient  variations  are 
shown  below.  These  are  presented  as  a  ratio  of  the 
scattering  coefficient  for  aerosols  to  the  scattering 
coefficients  for  aerosols  at  20 7.  relative  humidity  as  a 
function  of  relative  humidity. 

4 

3 

.MMH.aoo  2 

1 
o 

RH  in  X 

Figure  2-6  Scattering  Coefficient  Variation 


Models  of  the  aerosol  distribution  in  the  marine 
environment  can  be  very  complex  because  of  the  amount  and 
size  of  the  particles  present  in  the  marine  environment. 
These  are  functions  of  particle  generation  at  the  sea 
surface  and  the  vertical  transport  of  the  particles  once 
they  are  suspended  in  the  atmosphere. 

Wind  speed  dependence  on  aerosol  distribution  in  the 
maritime  model  is  similar  to  that  of  continental  models 
until  the  wind  speed  nears  14  knots.  At  14  knots  the 
aerosol  concentration  has  increased  by  a  factor  of  about  1.5 
due  to  aerosol  generation  by  the  surface  whitecaps  which 
form  at  this  wind  speed.  As  the  wind  speed  increases,  the 
concentration  may  rise  to  a  value  of  10  to  100  times  the 
calm  atmospheric  value  near  the  surface,  but  decrease  in 
concentration  with  altitude  up  to  approximately  5 
ki  1  ometers. 

It  can  now  be  seen  that  the  atmospheric  conditions  have 
a  dramatic  effect  upon  the  extinction  and  absorption  at 
various  wavelengths.  In  this  paper  the  interest  lies 
primarily  in  the  infrared  region  since  this  is  where  our 
model  (SIREM)  is  designed  to  simulate  target  acquisitions. 
Hence  it  is  good  practice  to  become  familiar  with  those 
factors  which  affect  the  systems  being  modeled. 

Incorporation  of  "real -world"  occurrences  is  a  vital  part  of 
systems  modeling.  Those  models  which  use  random  (non- 
deterministic)  methods  to  emulate  a  patterned  parameter  can 
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only  be  termed  less  than  valid.  They  Mould  serve  useful  for. 
nothing  more  than  gaming  exercises  since  the  answers  could 
not  be  relied  upon. 

Since  the  SIREM  model  is  designed  to  obtain  measures  of 
effectiveness  of  IR  decoy  deployment  versus  an  IR  seeking 
missile,  those  measures  may  only  be  obtained  through 
realistic  incorporation  of  the  factors  that  affect  the 
engagement.  The  primary  intent  of  the  folloMing  chapters  is 
to  define  the  scenario,  present  the  principles  of  detection, 
and  show  how  the  LOWTRAN  program  can  be  used  to  develop  a 
very  compact  and  accurate  method  of  incorporating 
atmospheric  absorption  into  the  acquisition  process. 


As  with  any  detection  method,  whether  RF,  visible,  or 
infrared,  it  is  assumed  that  the  target  is  either  emitting 
or  reflecting  energy  at  some  power  density  in  the  frequency 
band  of  the  detector  as  in  the  scenario  shown  below. 

i 


From  extended  ranges  this  emission  or  reflection  can  be 
assumed  to  appear  as  a  .point  source,  however  at  closer 
ranges  the  physical  size  of  the  target  must  be  expanded  into 
a  distributed  source  emitter,  especially  considering  the 
typical  optical  magnification  performed  by  the  seeker. 

Since  this  document  deals  specifically  with  infrared 
emission  and  transmission,  the  case  of  target  reflection 
will  henceforth  be  disregarded. 

The  unit  of  measure  given  to  the  power  emitted  from  the 
target  is  the  radiant  intensity  (J) ,  commonly  stated  in 
watts  per  steradian  (solid  angle).  The  target,  assumed  to 
be  a  point  source,  emits  a  source  level  radiation  in 
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specific  spectral  bands.  The  radiant  flux  density  of  any 
wavelength  across  the  continuum  has  a  natural  reduction  in 
value  based  on  the  inverse  square  law  with  range  from  the 
ship  (target)  to  the  missile  (receiving  platform). 

The  detector  in  the  seeker  head  is  composed  of  a 
material  which  electrically  (via  photon  energy)  sensitive  to 
incoming  radiation  across  some  defined  frequency  bandwidth 
inherent  to  the  composition.  This  bandwidth  is  preferrably 
fairly  narrow.  If  it  is  allowed  to  become  excessively  wide, 
the  signal-to-noise  ratio  is  reduced  an  unacceptable  level 
of  false  alarm  detections  as  discussed  by  Hudson  [Ref.  2:ch. 
123.  For  this  reason  most  detector  systems  have  external 
spectral  filters  to  narrow  the  bandwidth  and  eliminate  these 
mi sdetections. 

Based  on  the  composition  of  the  detector  material  each 
detector  will  have  an  optimum  electrically  active  wavelength 
to  which  it  is  sensitive.  As  the  incoming  energy  shifts  in 
frequency  to  either  side  of  the  optimum,  the  electrical 
activity  in  the  detector  is  reduced,  creating  some 
distribution  which  reflects  the  electrical  activity,  or 
molecular  transition  density  as  a  function  of  wavelength. 
This  relationship  should  be  known  in  order  to  reflect  any 
spikes  in  the  sensitivity  curve  where  certain  energy  bands 
may  be  particularly  absorptive  of  incident  energy.  The 


profile  shown  in  Figure  3-2  is  one  which  could  be  used  in 
the  typical  simulation  environment  based  on  the  above 
principles. 


Detector 
Radi  at  i  on 
Detectivity 
(Relative) 

2.0  3.0  4.0  X.. 

Figure  3-2. Detect ivi ty  over  Bandwidth 

The  actual  electrical  bandwidth  of  a  detector  has  three 
primary  dependencies:  first*  the  material  from  which  it  is 
constructed;  second,  the  temperature  at  which  it  is 
operated;  third,  the  physical  area  of  the  detector  surface. 
The  detectivity  (D~)  is  the  convenient  measure  categorizing 
a  detector's  operation  within  these  parameters.  It 
normalizes  the  "sensitivity"  of  various  detectors  to  an 
electrical  bandwidth  of  1  Hz,  and  an  area  of  1  square 
centimeter,  and  is  a  convenient  means  for  comparing  the 
equivalent  electrical  response  between  detectors  in  circuits 
of  different  bandwidth.  Shown  in  Figure  3-3  is  *  typical 
Lead  Sulfide  detector's  D~  character i st ics  taken  from  Hudson 
CRef.  2:pp.  3651. 
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Figure  3-3.  Detectivity  versus  Wavelength 


The  value  o-f  D*  is  dependent  upon  the  temperature, 
wavelength  and  bias  current  applied.  The  low  temperatures 
are  required,  as  well  as  bias  current,  to  maximize  the 
system's  si gnal -to-noi se  ratio. 

The  object  o-f  our  simulation  is  to  characterize  the 
search  and  acquisition  process  as  closely  to  actual  as 
possible.  To  do  this  we  must  model  the  detection  process 
(seeker),  atmospheric  effects,  flight  character i sties  of  the 
missile,  and  target  (ship  or  decoy)  signatures  as  closely  as 
possible.  Taking  all  factors  into  account,  the  maximum 
range  of  the  detec  tear  must  be  calculated. 

The  spectral  irradiance  from  a  target  (assumed  distant 
in  the  field  of  view  (FOV)  is 

Hv  =»  Jv  T.(V)  in  Watts/cm*  (Eq.  3.1) 

R* 


where  T»(\>  is  the  transmittance  o-f  the  path,  and  the 
subject  o-f  -following  chapters.  R  is  the  distance  from  the 
seeker  to  the  target. 

The  spectral  radiant  power  into  the  detector  is 

Px  =  HxA.T.(X)  in  Watts  (Eq.  3.2) 

where  Ao  is  the  area  of  the  entrance  aperature  of  the 
optics,  and  Ta<\)  is  the  transmittance  of  the  whole  sensor 
system.  T.  <\>  is  the  product  of  the  transfer  functions  for 
the  detector  and  all  external  filter  functions  such  as 
protective  windows,  lenses,  reticles,  and  maskings. 

The  signal  voltage  from  the  detector  is 

V.  =  Pxft<X>  (Eq.  3.3) 

where  H<X)  is  the  spectral  responsivity  of  the  detector. 

Thus  far  the  equations  represent  an  infinitesimal 
spectral  interval  centered  about  some  wavelength  X-  Since 
there  is  a  finite  bandwidth  to  an  actual  system,  this  should 
be  taken  into  account  as  shown  in  Hudson  CRef.  2:ch.  133. 
Remember,  however,  that  we  are  assuming  "worst  case"  in  the 
simulation  process.  So,  by  selecting  the  optimum  detectivity 
wavelength  of  the  particular  detector,  we  have  effectively 
maximized  the  range  of  detection  for  our  seeker.  This 
assumption  may  only  be  made  if  the  seeker  system  being 
modelled  has  a  relatively  narrow  spectral  response  due  to 
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external  -filtering  -for  a  better  si  gnal -to-noi  se  ratio.  This, 
is  virtually  always  the  case  -from  design  inception  since  the 
false  alarm  problem  is  critical  here. 

Further  development  of  the  equations  above  is  performed 
by  Hudson  and  leads  to  the  Optical  Trade-Off  Equation: 


R  =  [  JT*3‘/2 

!  f  D.  (NA)  T-l 1,2  [D~3*'2 

1  1 1/2 

1  2  1 

(*Af ) 1/2  ( V*/V«)  j 

\_1_/  \ _ 2 _ /  \_3_/  \ _ 4 _ /  (Eq.  3.4) 

1  -  Target  irradiance  and  atmospheric  transmittance 

2  —  Optics  parameters 

3  —  Detector  character i st ics 

4  -  System  parameters  and  signal  processing 

where  R  =  detection  range. 

D„  =  diameter  of  entrance  of  aperture  optics. 

NA  =  numerical  aperture  [Ref.  2:ch.  53. 

W  =  instantaneous  field  of  view  of  sensor. 

Af  =  frequency  bandwidth  in  Hz. 

V./Vm  =  signal  to  noise  ratio  [Ref.  2:pp.  4193. 

By  breaking  the  equation  into  sections  as  shown  above  we 
see  in  the  first  term  that  reduction  of  target  signature 
<ship)  and  enhancement  of  the  decoy  signature  inband  would 
shift  the  acquisition  probability  more  toward  the  decoy.  In 
the  engagement  process  nothing  can  be  done  about  the 
atmospheric  transmittance,  although  it  does  have  dramatic 
effects  on  the  acquisition  range.  Knowing  the  parameters  of 
the  second,  third  and  fourth  term  facilitates  modelling 
various  missile/seeker  systems. 


At  ranges  beyond  the  sensitivity  of  the  seeker  head, 
i.e.  where  incident  energy  is  insufficient,  detection 
probability  would  be  zero.  In  addition  to  the  inverse 
square  law  power  loss,  the  atmospheric  losses  previously 
discussed  must  be  taken  into  account  in  order  to  accurately 
determine  detection  probability.  Since  the  detectivity  of 
the  seeker  is  independent  of  the  atmospheric  effects  on  the 
emitted  signature,  the  resultant  probability  of  detection 
would  be  the  product  of  probability  of  detection  assuming  no 
atmospheric  effects,  and  the  atmospheric  transmittance  which 
varies  from  0  to  1. 

Few  models  if  any  actually  emulate  the  signal 
integrating  principle  upon  which  most  detectors  depend  for 
threshold  detection.  The  scanning  technique,  whether  serial 
or  parallel,  performs  the  temporal  integral  of  the  incoming 
energy  across  the  "window  function"  which  defines  the 
acceptance  bandwidth  of  the  seeker.  The  seeker  system's 
probability  of  detection  cur ve-versus— wavelength  is 
proportional  to  the  detector's  detect i vi ty— versus-wavelength 
curve. 

The  present  version  of  the  SIREM  model  emulates  a 
scanner  type  seeker  with  limited  field  of  view  typical  of 
the  ASMs  previously  discussed.  This  model  is  presently 
being  upgraded  to  account  for  atmospheric  losses  over  range. 
It  will  do  so  by  assuming  the  "worst  case"  posture.  By 
this,  it  is  meant  that  at  the  optimum  detector  wavelength. 
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the  model  presently  calculates  probability  of  detection  for 
incremental  missile  movements.  This  value  is  calculated  as 


a  "worst  case"  value  for  the  ship  based  on  the  detectivity 
of  the  seeker  and  the  maneuvering  capability  of  the  missile. 
By  multiplying  the  probability  of  detection  (assuming  no 
losses)  at  each  range  increment  by  the  corresponding 
transmittance  value  at  that  range,  the  target  acquisition 
probability  is  more  accurately  represented.  This  is  the 
case  whether  the  acquired  target  is  a  decoy  or  the  ship. 

As  an  example  assume  that  the  transmittance  has  been 
calculated  for  a  given  atmospheric  profile  for  the 
engagement  scenario.  The  transmi ttance  is  then  known  at 
each  increment  of  range.  As  the  missile  approaches  its 
target  the  atmospheric  transmittance  increases  approximately 
exponentially.  If  the  probability  of  detection  is 
calculated  for  the  seeker  system  in  each  incremental 
movement  performed  by  the  simulation,  the  engagement 
scenario  probability  of  detection  for  the  missile  would  be 
the  product  of  transmittance  and  original  seeker  acquisition 
probability  generated  by  its  own  separate  technique.  Hence 
if  the  original  P„  were  0.8,  and  the  transmittance  were  0.6, 
the  total  system  probability  of  detection  would  then  be 
0.48. 
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IV.  ATMOSPHERIC  MODELLING 


The  previous  chapters  presented  the  principles  of 
propogation  of  infrared  radiation  throughout  the  atmosphere 
and  the  factors  which  affect  the  acquisition  process.  This 
chapter  will  deal  with  an  empirical  approach  to  modelling 
the  atmospheric  transmi ttance  using  very  accurate  single 
equations  once  the  conditions  of  visibility,  range  and  wind 
speed  are  given. 

Atmospheric  transmittance  losses  for  infrared  radiation, 
not  unlike  radio  frequency  transmission  losses,  exhibit  the 
characteristic  1/ (Range)2  loss  factor  as  shown  below. 


1.0 

Atmospheric 
Transmi ttance 


0  5  10  15  20  Range  (Km) 

Figure  4-1.  Transmittance  versus  Range. 

An  actual  curve  varies  widely  in  its  rate  of  transmission 
loss  with  range  depending  heavily  upon  the  particular 
infrared  wavelength  in  question,  the  amount  and  size  of 
particulate  matter  in  the  air,  the  altitude(s),  and  weather 
conditions  as  discussed  previously.  The  LOWTRAN  program 
CRef.  II  was  developed  to  a  large  extent  from  empirical 


observations  and  is  used  to  predict  atmospheric 
transmittance  for  any  given  set  of  conditions.  The  desired 
goal  in  the  continuing  development  of  the  SIREM  model  was  to 
find  a  simple  technique  for  incorporating  the  calculation  of 
atmospheric  transmittance  into  the  incremental  range 
movements  of  a  simulated  incoming  missile  and  its  associated 
infrared  seeker  head  in  order  to  more  accurately  portray  the 
probability  of  acquisition  of  the  seeker. 

A.  LOWTRAN 

The  LOWTRAN  code  calculates  both  atmospheric 
transmittance  and  radiance  throughout  the  infrared  spectrum 
(0.25  to  28.5  micrometers) .  It  uses  a  single  parameter  band 
model  for  molecular  absorption,  and  includes  the  effects  of 
continuum  absorption,  molecular  scattering,  and  aerosol 
absorption.  Refraction  and  earth  curvature  are  included  in 
the  calculation  for  slant  paths. 

Five  seasonal  models  and  the  1962  U.S.  Standard 
Atmosphere  are  provided  as  inputs  to  the  LOWTRAN  program 
using  typical  altitude,  pressure,  temperature,  water  vapor 
density  and  ozone  density  profiles  for  each.  The  five  other 
atmospheric  profiles  consist  of  (1)  Tropical  (IS9  N)  ,  (2) 
Midlatitude  Summer  (45°  N-July),  (3)  Midlatitude  Winter  (45° 
N-January) ,  (4)  Subarctic  Summer  (60*  N-July>,  and  (5) 

Subarctic  Winter  (60°  N-January >.  With  these  profiles  the 
user  may  select  the  general  atmospheric  conditions 


applicable  to  the  particular  simulation.  Additionally  the 
user  may  input  radiosonde  data  if  desired  to  make 
calculations  for  a  specific  profile.  Particular  information 
regarding  mixing  ratios  and  assumptions  made  by  the  LOWTRAN 
program  can  be  obtained  from  the  LOWTRAN  Manual  [Ref.  1]  and 
the  associated  texts  listed  therein. 

The  LOWTRAN  code  includes  three  boundary  layer  aerosol 
types:  Rural,  Urban,  and  Maritime.  Due  to  the  nature  of  the 
SIREM/ASM  model,  only  the  latter  Mas  used  since  the 
composition  and  distribution  of  aerosols  of  oceanic  origin 
are  peculiar  to  the  IR/ASM  problem.  These  aerosols  are 
largely  sea— salt  particles  in  the  lower  boundary  layer 
caused  by  evaporation  of  sea— spray  droplets.  Together  with 
a  background  aerosol  of  more  or  less  pronounced  continental 
characteristics  they  form  a  fairly  uniform  maritime  aerosol 
which  is  representati ve  of  the  lower  2  to  3  kilometers  of 
the  atmosphere  over  the  oceans. 

These  aerosols  should  not  be  confused  with  the  heavy 
direct  sea-spray  aerosols  found  in  the  lower  10  to  20  meters 
above  the  ocean  surface  which  is  strongly  dependent  on 
recent  wind  velocities  and  sea  surface  turbulence.  The 
extinction  and  absorption  coefficients  are  calculated  as  a 
function  of  the  relative  humidity  and  based  heavily  upon  the 
refractive  index  for  the  given  altitude  and  weather 
conditions. 
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The  ability  to  calculate  slant-path  transmi ttances  is 
incorporated  utilizing  a  data  base  of  seasonal  and  latitude 
dependencies  of  vertical  aerosol  distributions.  This  allows 
transmittance  calculations  of  “high-flyer"  missile  attack 
profiles  or  aircraft  infrared  missile  engagements  by 
changing  the  input  parameters  to  the  LOWTRAN  program  to 
match  those  of  the  simulation  requi resents.  Iterating  this 
calculation  throughout  the  slant  path  range  values  yields 
the  data  points  necessary  to  model  the  reduction  in 
transmittance  over  the  entire  engagement— range  scenario  with 
a  single  equation.  The  horizontal  path  used  by  the  SIREN 
model  for  simulating  cruise  missiles  is  assumed  to  be  a 
constant  pressure  path  where  earth  curvature  and  refraction 
effects  are  negligible,  or  taken  into  account  through  the 
use  of  the  4/3  earth  radius  model .  An  alternative  flat 
earth  may  be  assumed  for  short  range  situations. 

The  LOWTRAN  code  calculates  transmittance  as  a  function 
of  "absorber  density"  for  the  path,  the  pressure,  and 
temperature,  for  the  particular  wavelength  band  chosen.  It 
utilizes  both  empirical  laboratory  data  and  available 
molecular  line  constant  data  in  performing  its  calculations. 
The  output  for  a  typical  LOWTRAN  transmittance  run  is  shown 
in  Figure  4-2.  The  input  to  the  program  for  spectral  range 
is  in  the  form  of  wavenumbers  (inverse  centimeters).  In 
this  example  the  range  is  900  to  1145  cm*1  at  increments  of  5 
cm*1  required  by  the  LOWTRAN  program.  The  lesser  number  of 
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Figure  4.2.  LCWTHAN  Program  Output 


wavenumbers  is  always  the  -first  input  into  the  LOWTRAN 
program,  while  the  higher  wavenumber  is  the  second  input. 
This  corresponds  to  the  high-wavelength  to  low-wavelength 
output  in  column  2  o-f  Figure  4-2.  The  conversion  to 
wavenumbers  may  be  made  using  the  -following  -formulas. 

X  =  lO"*/  (Wavenumbers)  or 

(Eq.  4.1) 

Wavenumbers  =  10'VX 

Here  the  wavenumbers  are  in  inverse  centimeter  values  as 
used  by  the  LOUTRAN  program.  Transmittance  values  -for  each 
o-f  the  contributing  absorbers  (H20,  C02,  Ozone,  N2,  Molecular 
Scattering,  Aerosol,  and  Nitric  Acid)  are  then  calculated 
-for  each  wavenumber  (wavelength  in  column  2)  in  row  format. 
The  total  transmittance  for  each  wavenumber  is  them 
calculated  as  the  product  of  the  transmi ttances  of  each  of 
the  absorbing  elements  and  is  listed  in  column  3.  The 
average  transmittance  listed  at  the  bottom  of  Figure  4-2  is 
simply  the  average  of  the  transmi ttances  in  column  3,  and 
reflects  the  "integrated"  overall  transmittance  across,  in 
this  case,  the  8.7  to  11.1  micrometer  bands  for  the  chosen 
input  parameters  (model,  visibility,  range,  etc.). 

Since  the  bandwidth  of  most  IR  detectors  is  limited, 
either  by  its  own  material  limitations  or  external  filters, 
the  input  wavenumber  range  will  seldom  be  as  wide  as  the 
previous  example.  The  3.5  micrometer  band  was  selected  for 
the  development  of  this  paper  and  consists  of  wavenumber 
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inputs  o-f  2850  (lower)  to  2860  (upper)  cm*1,  and  correlates 
to  3.496-3.509  micrometers  in  5  cm*1  intervals. 

In  addition  to  the  output  shown  in  Figure  4-2  the 
program  generates  a  large  quantity  o-f  "header**  information 
for  each  run,  which  reprints  the  inputs  to  the  program.  The 
unfamiliar  user  should  obtain  hardcopy  of  early  runs  to  see 
what  the  output  consists  of.  Once  comfortable  with  the 
location  of  various  values  of  the  printout,  use  of  the  video 
screen  is  more  appropriate  and  less  wasteful  of  paper.  The 
objective  of  each  run  is  the  “average  transmittance"  for  the 
given  input  parameters,  which  will  later  be  used  in  the 
curve  fitting  process. 

One  can  easily  see  that  in  order  to  curve  fit  the 
transmittance  as  a  function  of  range  with  each  other 
parameter  held  constant,  enough  runs  through  the  LOWTRAN 
program  must  be  conducted  at  "selected  ranges"  to  obtain  a 
smooth  curve  fit.  This  means  that  for  every  model,  range, 
wind  condition,  etc.  the  LOWTRAN  program  must  be  run  in 
order  to  obtain  the  appropriate  transmi ttance  value.  This 
would  be  a  monumental  task  for  the  computer,  and  without 
stripping  the  PRINT  statements  from  the  LOWTRAN  code,  could 
generate  a  massive  output  file.  The  objective  here  is  to 
show  that  a  minimum  number  of  selected  ranges  (7  or  8)  can 
be  used  for  this  process  to  decrease  computer  time  spent. 

The  task  becomes  much  easier  if  a  particular  set  of 
conditions  is  in  mind,  i.e.  a  Midlatitude  Summer  model  with 


average  wind  conditions  of  5  m/s  and  visibility  of  30 


kilometers  for  example.  This  circumstance  could  be  modelled 
by  running  a  "stacked  input"  facility  of  the  (-OUTRAN  program 
containing  the  ranges  desired.  A  test  run  can  be  performed 
to  find  the  approximate  maximum  range  where  the 
transmi ttance  value  falls  below  some  operator  desired 
threshold;  e.g.  T=0.0025. 

Curve  fitting  atmospheric  transmittance  as  a  function  of 
range  requires  that  all  other  parameters  (i.e.  wind, 
visibility,  etc.)  remain  constant.  Additionally,  enough 
points  on  the  range  scale  must  be  taken  to  provide  the 
necessary  accuracy  across  the  span  of  range  of  interest. 

The  missile  engagement  scenario  requires  that  transmittance 
be  a  factor  in  the  acquisition  process  from  horizon  to 
target,  thus  the  span  in  range  should  be  from  approximately 
20  kilometers  to  zero.  After  extensive  empirical  analysis 
the  values  in  Table  4.1  were  selected  to  meet  the  above 
cri teria. 

Up  to  this  point  it  was  assumed  that  all  parameters 
(including  wind  and  visibility)  were  constant.  The  desired 
goal  is  to  analyze  the  individual  effects  of  wind  and 
visibility  on  transmittance  and  develop  an  empirical 
equation  which  generates  a  good  approximation  to  the 
transmittance  value  for  any  wind  and  visibility  condition 
over  the  scenario  range,  similar  to  that  shown  in  Figure  4— 
1.  The  transmi ttance  takes  on  a  negative  exponential 
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characteri stic  whose  decay  constant  varies  as  a  -function  of 
wind  and  visibility.  The  idea  is  to  empirically  establish 
the  decay  constant  as  a  function  of  both  wind  and 
visibility,  the  two  principal  parameters  which  must  remain 
variable  to  the  model  user  to  real istical ly  model  any 
engagement  scenario.  Intuitively  (and  verified  by  LOWTRAN) 
atmospheric  transmittance  decays  much  more  rapidly  if  either 
visibility  decreases  or  wind  increases.  To  establish  how 
the  decay  constant  changes  with  wind  and/or  visibility, 
enough  values  for  both  wind  and  visibility  must  be  analyzed 
to  curve  fit  the  decay  constant  dependence  to  either.  The 
values  of  wind  and  visibility  in  Table  4.1  were  empirically 
selected  and  adequately  cover  the  conditions  for  most 
scenarios.  Then,  for  each  combination  of  wind  and 
visibility,  transmittance  values  over  the  range  span  are 
obtained  using  LQUTRAN.  This  builds  a  3— dimensional  table 
which  is  tha  basis  for  curve  fitting  and  establishing  the 
decay  constant  as  a  function  of  three  variables;  i.e.  range, 
visibility,  and  wind. 

TABLE  4.1.  RECOMMENDED  CURVE  FIT  DATA  POINTS 
Range  i n  Km  Mind  m/s  Visibility  in  Km 


0.2 

0.  1 

40 

1.0 

5.0 

30 

2.0 

10.0 

20 

4.0 

15.0 

15 

8.0 

20.0 

10 

12.0 

18.0 

21.0 

r 

f 
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The  method  used  -For  data  reduction  is  an  extremely  simple 
technique  implemented  For  microcomputer  spreadsheets.  The 
method  employs  a  scheme  For  perForming  Weighted  Least 
Squares  curve  Fit  techniques  to  the  data,  and  while 
calculating  the  coeFFicients  oF  the  curve  Fit,  the  goodness 
oF  Fit  oF  the  curve  to  the  data  may  be  viewed  graphically  by 
the  modeller  For  approval  or  correction  as  necessary.  Some 
examples  oF  the  process  will  be  given  later. 

The  LOWTRAN  program  serves  as  an  excellent  tool  For 
calculating  transmittance  over  a  wide  variety  oF  parameters, 
however,  the  program  consists  oF  9000  lines  oF  FORTRAN  code. 
Most  simulations  are  already  quite  large,  and  the  addition 
oF  this  code  size  as  a  subroutine  is  undesireable, 
considering  that  each  increment  in  range  oF  the  missile 
Flight  path  must  vector  to  LOWTRAN  as  a  subroutine  in  order 
to  obtain  an  updated  single  transmi ttance  value. 

The  alternative  method  presented  here  For  obtaining  the 
transmittance  value  desired  is  by  utilizing  LOWTRAN  external 
to  the  simulation  and  curve  Fitting  the  output  data  points 
as  a  Function  oF  range  to  a  single  equation,  given  the 
visibility  and  wind  conditions  oF  the  particular  engagement. 
This  method  requires  extensive  use  oF  the  LOWTRAN  program  to 
obtain  the  data  points  necessary  to  perForm  the  curve  Fit. 
Once  completed  though,  the  two  coeFFicients  oF  the 


approximated  curve  are  the  only  data  necessary  to  write  an 
equation  of  the  form: 

T  =  0  Exp(-«R)  <Eq.  4.2) 

The  coefficients  of  this  equation  (0  and  a)  are  all  that  is 
necessary  to  allow  calculation  of  transmi ttance  as  the 
missile  moves  incrementally  through  the  stochastic 
simulation.  When  an  increment  inrange  occurs,  a  vector  to 
the  transmittance  calculating  subroutine  (or  function)  with 
this  equation  updates  the  transmittance  affect  on  the 
probability  of  acquisition.  As  previously  discussed,  the 
decay  constant  (a)  is  actually  a  function  of  wind  and 
visibility.  Later  in  this  chapter  accomodations  will  be 
made  for  this. 

The  SIREM  model  specifically  calculates  through 
simulation,  the  effectiveness  of  deployment  of  Torch/PIP 
Torch  type  infrared  decoys  against  infrared  seeking  anti¬ 
ship  missiles  (ASMs) .  It  is  assumed  from  this  point  on  that 
the  incoming  missiles  utilize  a  "low  altitude”  attack 
profile  typical  of  the  cruise  missile  variety  which  prevents 
early  detection  by  the  target  ship.  Assume  for  the  time 
being  that  it  flies  at  an  altitude  of  approximately  60  feet 
(20  meters  from  the  LOWTRAN  example  modelled  herein).  This 
would  constitute  a  horizontal  transmi ttance  path  when  using 
the  LOWTRAN  program  to  calculate  transmi ttance  throughout 
the  missi le-to-target  flight  path.  Some  other  assumptions 
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utilized  in  the  empirical  derivations  in  this  paper  which 


were  used  as  inputs  to  the  LOWTRAN  program  are: 

a)  The  Marine  Atmospheric  Boundary  Layer  model  was  the 
only  one  utilized. 

b)  A  land  haze  -factor  o-f  3  was  used  (default). 

c)  The  wavelength  modelled  here  was  3.5  micrometers. 

d)  Each  of  the  6  atmospheric  profiles  was  calculated. 

e)  Default  parameters  inherent  to  LQWTRAN  were  used, 
varying  only  the  visibility  conditions,  range,  and 
wind  conditions 

f)  There  was  no  precipitation. 

Each  of  these  parameters  may  be  changed  if  desired  to 
meet  an  alternative  modelling  requirement.  For  example,  if 
it  is  known  that  a  “new"  IR  missile  threat  has  evolved  which 
utilizes  a  "high-flyer"  attack  profile  and  a  3.2  micrometer 
optimized  detector,  the  modeller  could  then  choose  the 
slant-path  option  in  the  LOWTRAN  program  selected  at  the 
necessary  wavelength.  Knowing  only  the  flight  profile, 
detector  f requency ( i es)  and  atmospheric  conditions  would 
then  allow  full  use  of  the  empirical  transmittance  modelling 
technique  presented  here  without  insertion  of  9000  lines  of 
LOUTRAN  code. 

Another  advantage  to  this  approach  other  than  its 
simplicity  and  accuracy  is  that,  if  the  simulation  is 
written  in  a  language  other  than  FORTRAN,  this  method  can 
still  be  easily  adapted  due  to  the  short  code  length 


required  using  the  standard  exponential  -function.  The 


equational  -form  is  easily  translated  as  will  be  seen. 

Let's  analyze  an  example  scenario.  Suppose  we  are  faced 
with  our  ship-versus-missi le  engagement  described 
previously,  and  wish  to  calculate  atmospheric  transmittance 
at  various  ranges  in  order  to  approximate  its  affect  on  the 
acquisition  process.  First  we  must  assume  some 
meteorological  conditions  appropriate  to  our  example  for 
inputs  to  the  LOWTRAN  program.  We  will  select  here  the 
parameters: 

1)  Mid latitude  Summer  profile. 

2)  Visibility  is  30  kilometers. 

3)  A  maritime  environment. 

4)  A  land  haze  factor  of  3  (may  vary  from  1  to  10,  with 
lO  being  bad  haze  as  around  Los  Angeles). 

5)  Wind  is  virtually  nill.  (0.1  m/s)  and  has  been  so 
for  the  past  24  hours. 

6)  The  range  from  the  seeker  to  target  is  4  kilometers. 

7)  Assume  a  3.5  micrometer  band  (2850-2860  cm-1) 

8)  The  path  is  horizontal  at  20  meters  altitude. 

We  may  now  insert  out  data  into  the  LOWTRAN  program  in  the 
format  listed  in  the  LOWTRAN  Manual  [Ref.  11. 

In  doing  so  we  would  find  that  at  our  4  kilometer  range, 
the  transmittance  for  this  wavelength  would  be  0.6516  (See 
Appendix  A:  Table  A. 2. a).  The  value  may  differ  only 

This  means  that  the  received 
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slightly  between  computers 


energy  density  at  a  range  of  4  kilometers  is  only  65%  as 


high  as  that  for  a  perfectly  transmitting  medium  (vacuum). 
We  can  now  duplicate  our  efforts  for  the  4  kilometer  range, 
but  now  with  a  wind  velocity  of  10  meters  per  second, 
keeping  visibility  the  same.  Doing  this  we  obtain  a 
transmittance  value  of  0.4537.  In  the  diagram  below  we  see 
that  increasing  the  wind,  or  decreasing  visibility,  has  a 
dramatic  effect  on  the  decay  constant  associated  with  the 
transmi ttance. 


Atmospheri c 
Transmi ttance 


Range  (Km) 


Figure  4-2.  Wind  Effects  on  Transmittance. 


Note  the  approximately  exponential  fall  off  with  range. 
This  is  the  curve  we  wish  to  approximate  as: 

T(R,V,W)  =  tf -EXP (« ( V) *f (W) *R • 6 )  (Eq.  4.3) 

from  the  data  points  obtained  from  the  LOWTRAN  program. 
Here  T,  R,  V,  and  W  represent  transmi ttance,  range, 
visibility,  and  wind  respectively.  There  is  correlation 
between  wind  and  visibility,  and  the  effects  are  accounted 
for  in  the  factor  6. 


In  order  to  model  any  particular  atmospheric  profiles, 
the  data  points  for  the  3-dimensional  table  may  be  obtained 
by  adapting  the  following  overlay  algorithm  as  the  main 
program,  using  LOWTRAN  as  a  subroutine. 

FOR  Atmospheric_Prof i 1 e  -  1  to  6 
FOR  Wind_Speed  =  0  to  20  Step  5 
Wind  =  0.01  +  FLOAT <Wind_Speed) 

FOR  I  =  1  to  Nr_of _Range_Values 
READ  Range_Value 
OPEN  Input_File 
WRITE  Input_File 
CALL  LOWTRAN 
WRITE  Output 
NEXT  I 

NEXT  Visibility 
NEXT  Wind_Speed 
NEXT  Atmospher i c_Prof i 1 e 

Note  1:  Prior  to  entering  the  looping  structure  the  user 
must  assign  values  to  the  other  parameters  used  by  the 
LOWTRAN  program  (e.g.  wavenumbers,  haze  factor,  etc.). 

Note  2:  The  actual  wind  speed  may  not  be  zero.  If  so 
LOWTRAN  defaults  to  4.1  m/s.  Hence  the  following  statement. 
Note  3:  The  input  file  must  be  formatted  as  specified  in 


*  See  Note  1 

*  See  Note  2 


*  See  Note  3 


the  LOWTRAN  Manual  CRef.  11. 


Again  it  should  be  emphasized  that  the  source  code  o-f 
the  LOWTRAN  program  should  be  commented  out  to  prevent 
printing  o-f  the  massive  administrative  data.  Only  the 
“Average  Transmittance"  over  the  sufficiently  narrow  band  is 
desired.  The  output  may  be  directed  to  a  printer  or  file. 
The  latter  method  is  preferred  if  input  to  the  curve  fitting 
process  is  desired  in  a  direct  fashion.  In  order  to  keep 
track  of  which  data  value  is  associated  with  which  range, 
visibility,  and  wind  value,  it  may  be  wise  to  write  those 
associated  values  in  the  same  output  row. 

Continuing  with  our  example  we  will  assume  for  now  that 
only  one  value  for  wind,  and  likewise  visibility  are 
calculated.  The  transmittance  would  then  become  only  a 
function  of  range  since  all  other  parameters  remain 
constant.  The  desired  goal  is  to  select  range  values  which 
best  facilitate  the  curve-fitting  process.  Using  the  range 
values  from  Table  4.1  we  can  now  perform  an  example  of  the 
curve  fitting  process,  while  presenting  the  spreadsheet 
technique  used  for  doing  Weighted  Least  Squares. 


1 .  Ex  amp 1 e 


Making  the  assumption  that  all  parameters  remain 
constant  as  before,  with  only  range  changing,  a  curve  fit  of 
the  following  LOWTRAN  data  is  performed  for  the  Midlatitude 
Summer  profile  with  wind  at  5  meters  per  second,  and 
visibility  of  30  kilometers. 

TABLE  4.2a.  LOWTRAN  DATA  FOR  CURVE  FITTING 


RANGE 

TRANSMITTANCE 

WEIGHTING  FACTOR 

0.2 

0.9272 

8 

1.0 

0.7824 

1 

2.0 

0.6588 

2 

4.0 

0.4828 

21 

8.0 

0.2723 

6 

12.0 

O. 1588 

1 

18.0 

0.0727 

1 

21.0 

0. 0498 

2 

The  Weighted  Least  Squares  procedure  is  conducted  in 
two  parts.  To  effect  the  weighting,  implementing  multiple 
occurrence  of  the  same  data  points  is  shown  in  table  4.3. 
Following  in  figure  4.4a  and  4.4b  is  the  actual  Least 
Squares  algorithm  implemented  in  spreadsheet  format.  The 
output  desired  is  the  values  for  the  A-Coef f icient  (0)  and 
B-Coef f icient  <«)  for  equation  4.2.  In  this  example  the 
first  coefficient  4=0.875824,  and  the  second,  « =-0. 140899. 
Using  equation  4.2  yields  the  table  which  follows  at  the  top 
of  the  following  page. 
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TABLE  4.2b.  DATA  POINTS  EXTRACTED  FROM  FITTED  CURVE 


RANGE 

TRANSMITTANCE 

0.2 

0.8515 

1.0 

0.7607 

2.0 

0.6607 

4.0 

0.4985 

8.0 

0.2837 

12.0 

0. 1615 

18.0 

0.0693 

21.0 

0.0454 

Comparing  table  4.2a  with  4.2b  (i.e.  the  LOUITRAN 
data  versus  the  exponential  Curve-Fitted  data)  we  notice 
comparable  transmittance  values,  especially  at  ranges 
greater  than  one  kilometer.  Due  to  the  nature  of  the 
engagement  scenario,  we  may  assume  transmittance  values  at 
ranges  less  than  one  kilometer  or  so  to  be  equal  to  1.0, 
since  its  effect  on  acquisition  at  these  ranges  is  minimal. 
Above  1  kilometer  the  error  is  generally  one  to  two  percent 
at  most,  and  hence,  very  useable  for  "probability" 
modelling.  Within  the  following  tables  the  Least  Squares 
process  in  implemented  into  spreadsheet  format.  Although  in 
this  case  Lotus  Development  Corporation's  1-2-3  (TM)  was 
used,  the  technique  is  similar  if  using  commercially 
available  packages.  The  data  is  arranged  in  row— versus- 
column  format  with  the  Least  Squares  algorithm  implemented 
in  Tables  4.4a  and  4.4b. 


& 


A 


B 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


Model t 
Visibi 1 ity* 
Variad  paraa: 

Range  (km) 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

1 

2 

2 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

8 

8 

8 

8 

8 

8 

12 

18 

21 

21 


Midi at  Sue 
30 

Visibility 

Transmi ttance 
0.9272 
O. 9272 
0.9272 
0.9272 
0.9272 
O. 9272 
0.9272 
O. 9272 
O. 7824 
0.6588 
0.6588 
O. 4828 
0.4828 
0.4828 
0.4828 
0.4828 
O. 4828 
O. 4828 
O. 4828 
0.4828 
O. 4828 
O. 4828 
O. 4828 
O. 4828 
0.4828 
0.4828 
0.4828 
0. 4828 
0.4828 
0.4828 
0.4828 
0.4828 
0.2723 
0.2723 
O. 2723 
0.2723 
0.2723 
0.2723 
O. 1588 
0.0727 
0.0498 
0.0498 


TABLE  4.3. 

TRANSMITTANCE  AND  RANGE 
DATA  POINTS  USED 


NOTES 

Using  the  same  data  point 
multiple  times  at  any  point 
minimizes  error  about  that 
point  in  the  curve  fit 
process.  This  emulates  the 
weighting  matrix  commonly 
used  in  performing  Wieghted 
Least  Squares  curve  fitting. 


TABLE  4.4a.  CALCULATIONS  USED  IN  LEAST  SQUARED  FITTING 
:  C  s  D  s  E  i  F 


1: 

2: 

3i 

5 

sun  of  Rt 

210.6 

Ri* 

This  section  finds  the  coefficients 

A  and  B  for  the  equations 

T  *  <T*Exp(«*R) . at  Visibility  =30  km 

Ln(T,)  CLn(T,)3*  Rt*Ln(T,) 

6 

0.04 

-0.0756 

0.0057 

7 

0.04 

-0. 0756 

0.0057 

-0.0151 

8 

0.04 

-0.0756 

0.0057 

-0.0151 

9 

0.04 

-0. 0756 

0.0057 

-0.0151 

10 

0.04 

-0.0756 

0.0057 

-0.0151 

11 

0.04 

-0.0756 

0.0057 

-0.0151 

12 

0.04 

-0.0756 

0.0057 

-0.0151 

13 

0.04 

-0. 0756 

0.0057 

-0.0151 

14 

1 

-0.2454 

0.0602 

-0.2454 

IS 

4 

-0.4173 

O. 1742 

-0.8347 

16 

4 

-0.4173 

0. 1742 

-0.8347 

17 

16 

-0. 7282 

0.5302 

-2.9126 

18 

16 

-O. 7282 

0.5302 

-2.9126 

19 

16 

-O. 7282 

0.5302 

-2.9126 

20 

16 

-0. 7282 

0.5302 

-2.9126 

21 

16 

-O. 7282 

0.5302 

-2.9126 

22 

16 

-0. 7282 

0.5302 

-2.9126 

23 

16 

-0.7282 

0.5302 

-2.9126 

24 

16 

-0.7282 

0.5302 

-2.9126 

25 

16 

-0.7282 

0.5302 

-2.9126 

26 

16 

-0. 7282 

0. 5302 

-2.9126 

27 

16 

-O. 7282 

0.5302 

-2.9126 

28 

16 

-0. 7282 

0.5302 

-2.9126 

29 

16 

-0. 7282 

0.5302 

-2.9126 

30 

16 

-0.7282 

0.5302 

-2.9126 

31 

16 

-0. 7282 

0.5302 

-2.9126 

32 

16 

-0. 7282 

0.5302 

-2.9126 

33 

16 

-0. 7282 

0.5302 

-2.9126 

34 

16 

-0. 7282 

0.5302 

-2.9126 

35 

16 

-0. 7282 

0.5302 

-2.9126 

36 

16 

-0.7282 

0.5302 

-2.9126 

37 

16 

-0. 7282 

0.5302 

-2.9126 

38 

64 

-1.3009 

1 . 6922 

-10.4068 

39 

64 

-1.3009 

1 . 6922 

-10.4068 

40 

64 

-1.3009 

1 . 6922 

-10.4068 

41 

64 

-1.3009 

1 . 6922 

-10.4068 

42 

64 

-1.3009 

1 . 6922 

-10.4068 

43x 

64 

-1.3009 

1 . 6922 

-10.4068 

44 

144 

-1.8401 

3.3860 

-22.0813 

45s 

324 

-2.6214 

6.8718 

-47. 1855 

46 

441 

-2.9997 

8.9984 

-62.9945 

47 

441 

-2.9997 

8.9984 

-62.9945 

48 

Sum(Rt>s 

Sum  Ln  <  Ti ) 

Sum  Squares  Sum 

of  Products 

49  s 

2079.3200 

-35.2421 

49.9966 

-320.8972 

47 


E>! 


TABLE  4.4b.  SPREADSHEET  LEAST  SQUARES  CALCULATIONS 


G  s  H 


1: 

2: 

3: 

4: 

5: 

6: 

7: 

8: 

9: 


Scratchpad 
-O. 1409 
-0. 1326 
20.3153 
20.4251 

A-Coef f icient 
B— Coe-f  -f  icient 
R-Squared 


> 

> 

> 

> 

0.875824  > 
-0.140899  > 
0.9946  > 


Formula  in  cell  at  left: 
Note:  Here  n=42 

<42*F49-D49*C3> / (42*C49-C3""2) 
(D49— G2*C3) /42 
G2* (F49— C3*D49/42) 
E49-D49^2/42 

Exp (G3) 

G2 

G4/G5 


This  satisfies  the  Least  Squares  process.  Cells  G2  through 
G5  represent  respectively: 


Gl) 

EnE(R,ln  (T, 

)  ) 

-  El  n  (T, )  ER,  1  /  CnER,2 

G2) 

EEln  (T, )  - 

(Gl)  (ER, )  23/n 

G3) 

G2  CER,  1  n  (T, ) 

-  ER,Eln(T,)/nD 

G4) 

Eln(T,)*  - 

(El  n  (T,  )  )  */n 

( ER, )  2  3 


(Eq.  4.4—4. 7) 


Cell  C3  is  the  sum  of  the  range  values  of  column  "A". 

C3)  ER,  (Eq.  4.8) 

The  R-Squared  value  represents  the  "Goodness  of  Fit" 
of  the  function  to  the  original  data  when  using  the 
approximating  equation.  In  this  case  the  fit  is  a  very  one, 
but  in  cases  of  either  extremely  good  or  poor  visibility 
(i.e.  >40  or  <5  kilometers)  errors  between  LOWTRAN 
predictions  and  empirical  results  grow,  and  the  fitted  curve 
will  have  regions  where  error  may  be  on  the  order  of  five  to 
six  percent.  For  modelling  in  the  probabilistic  sense  this 
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is  still  adequate,  since  the  region  o-f  error  is  transient 
and  errors  of  this  magnitude  have  only  small  effects  on  the 


acquisition  probability.  Figure  4-3  is  a  plot  of  the  above 
example  using  both  LOWTRAN-obtained  data  points,  and  the 
corresponding  exponential  curve  fit  data  points  obtained  by 
using  the  spreadsheet  technique  presented.  Me  may  now 
extend  our  process  to  calculate  the  atmospheric 
transmittance  for  any  of  the  6  LOWTRAN  profiles  for  any  wind 
or  visibility  conditions. 


Figure  4  —  3 


Comparison  —  LOWTRAN  vs  Cxponsntlal  Fit 


Extending  the  Curve  Fit 


The  3-dimensional  table  (A. 2. a)  of  transmittance 
values  has  been  obtained  for,  in  this  case,  the  Midlatitude 
Summer  Profile.  Thus  far  we  have  completed  the  curve¬ 
fitting  process  for  one  particular  wind  and  visibility 
condition  (i.e.  visibility  of  30  Km,  wind  speed  5  m/s  as 
calculated  in  Table  4.4).  The  A  and  B  coefficients  (<f  and 
a),  which  were  derived  as  a  result  of  Least -Squares  fitting 
an  exponential  function,  were  used  in  equation  4.2  to 
emulate  extinction  over  range  for  this  wi nd/visibi 1 i ty 
combination  to  obtain  Figure  4-3.  There  the  accuracy  of  the 
exponential  fit  versus  the  original  LOWTRAN  data  is  visually 
apparent.  Tables  A.l.a  through  A.&.b  in  Appendix  A  list  the 
values  of  transmittance  obtained  from  the  LOWTRAN  program 
and  from  the  curve  fitting  process  for  each  of  the  six 
atmospheric  profiles  available.  The  curve-fit  coefficients 
for  the  other  five  atmospheric  profiles  are  contained  in 
Tables  A. 7  and  A. 8  in  Appendix  A. 

We  now  wish  to  extend  this  procedure  to  be  able 
to  emulate  extinction  over  range  for  any  wind/visibility 
combination.  To  do  so  we  must  first  build  a  table  of  the  A 
and  B  coefficients  and  examine  how  their  values  change  from 
one  wind/visibility  combination  to  another.  Using  the 
combinations  from  Table  4.1  the  coefficients  were  obtained 
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•for  the  entire  atmospheric  pro-file  o-f  our  example.  These 
are  contained  in  Table  4.5. 

By  examining  Table  4.5  me  see  that  the  A 
coefficients  hardly  change  throughout  entire  range  of  wind 
and  visibility  values.  This  coefficient  may  be  assigned  a 
constant  value  for  the  remainder  of  the  extinction  curve-fit 
process.  This  is  true  when  modelling  each  of  the 
atmospheric  profiles  individually;  however ,  the  A- 
coefficient  does  change  from  one  profile  to  another.  This 
can  be  seen  in  Table  A. 7  in  Appendix  A. 


TABLE  4.5.  EXPONENTIAL  CURVE  FIT  COEFFICIENTS 
FOR  "MIDLATITUDE  SUMMER"  PROFILE 


A-Coef f icient  Value 


Vi s=40  Vis=30  Vis=20  Vis=15  Vis=10 


Mind 

— 

O 

O.B75788 

0.875808 

0.875754 

0.875834 

0.875812 

5 

0.875747 

0.875824 

0.875699 

0.875477 

0.877011 

io 

0.875820 

0.875716 

0.875650 

0.874656 

0.874052 

15 

0.875753 

0.875685 

0.875366 

0.876011 

0.872420 

20 

0.875777 

0.875885 

0.875886 

0.877431 

0.859462 

B-Coef f icient  Value 

Mi  nd 

Vi s=40 

Vi s=30 

Vi s=20 

Vis=15 

Vi s=10 

0 

-0.064275 

-0.065945 

-0.069313 

-0.072711 

-0.079494 

5 

-0. 118579 

-0. 140899 

-O. 185840 

-0.231112 

-0.322254 

IO 

-0. 129858 

-0. 156395 

-0.209959 

-0. 263693 

-0.371488 

15 

-0. 132919 

-0. 160663 

-0.216479 

-0.273163 

-0.384626 

20 

-0. 134394 

-0. 162714 

-0.219772 

-0. 277832 

-0.386820 

The  B-coef f i ci ents  do  vary  throughout  Table  4.5. 
Through  further  examination,  the  change  in  its  value  is  non¬ 
linear  with  wind  and  visibility,  and  as  such,  simple  linear 


interpolation  of  this  table  is  not  accurate.  Nor  is  this 
technique  desirable  because  it  requires  many  more  B- 
coefficient  data  points  than  those  contained  in  Table  4.5  to 
be  stored  in  memory  to  be  even  moderately  accurate.  The 
process  of  extending  equation  4.2  to  the  more  general  form 
of  equation  4.3  becomes  a  matter  of  finding  out  how  the  B- 
coefficients  in  Table  4.5  change  as  a  function  of  wind  and 
visibility  (i.e.  find  the  exponential  decay  constant  of 
extinction  in  terms  of  its  dependence  on  wind  and 
visibility),  and  fit  functions  which  will  accomodate  the  2- 
dimensional  decay  constant  variation. 

Thus  the  curve  fit  process  is  continued,  only 
now  it  must  be  performed  on  the  B-coef  f  l  c  l  ents  which  vary 
with  both  wind  and  visibility.  This  will  leave  'Range'  as 
the  only  variable  to  input  to  equation  4.3  for  calculating 
transmi ttance.  From  table  4.5  it  can  be  seen  that  the  B- 
coefficient  increases  negatively  (decreases)  as  either  1) 
wind  increases,  or  2)  visibility  decreases,  and  is  therefore 
a  function  of  wind  and  visibility.  Equation  4.3  stipulates 
that  the  decay  constant  for  extinction  over  range  consists 
of  three  parts.  These  ares 

a)  «(V)  -  Due  to  changes  in  visibility, 

b)  f (W>  -  Due  to  changes  in  wind,  and 

c)  6  -  Any  correlation  between  the  first  two. 
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By  also  observing  that  the  B-coef f icient  is 
least  negative  where  wind  is  minimum  and  visibility  is 
maximum,  we  can  now  use  the  ratios  o-f  any  other  B- 
coeff icients  relative  to  the  first  and  find  their 
mul tipi icati ve  increase  over  the  original  (baseline  — 
maximum  visibility,  minimum  wind)  B-coeff icient  due  to  wind 
and  visibility.  This  done  by  dividing  each  value  in  the  B- 
coeff icient  table  by  the  'baseline'  value  in  the  upper  left 
corner  of  the  table  as  is  done  for  Table  4.6.  Table  4.6 
depicts  how  much  greater  (negatively)  the  decay  constant  is 
for  any  particular  wind  and  visibility  combination  than  it 
would  be  for  its  best  transmittance  (i.e.  wind  zero, 
visibility  maximum).  The  non-linearity  of  the  B-coef f i ci ent 
with  wind  or  visibility  is  much  more  evident  in  Table  4.6. 

We  shall  see  in  a  few  moments  why  this  table  is  necessary. 


TABLE  4.6.  RATIO  OF  B(Vis,Wind)  TO  B  ( Vi  s„, ,  Wind*,*) 


Wi  nd 

0 

Vis 

40 

Vis 

30 

Vis 

20 

Vis 

15 

Vis 

10 

1.0 

1 . 02598 

1 . 07838 

1. 13125 

1 . 23678 

5 

1 . 84487 

2. 19213 

2.89133 

3.59567 

5.01367 

10 

2.02035 

2.43322 

3.26657 

4. 10257 

5.77966 

15 

2.06797 

2.49962 

3.36801 

4.24991 

5.98407 

20 

2.09092 

2.53153 

3.41924 

4.32255 

6.01820 

The  following  steps  are  the  procedures  used  to 
obtain  the  three  functions  of  the  exponential  decay  constant 
of  equation  4.3.  Individually  these  steps  produce  the 


functions  « (V) ,  f (W) ,  and  <5  which  will  allow  the  input  of 


any  wind  or  visibility  value  -from  0-20  m/s  or  40-5  Km 
respectively,  and  generate  the  approximated  B-coef f icient  of 
extinction  for  that  condition.  Keep  in  mind  that  this 
example  only  covers  the  ‘Midlatitude  Summer'  profile.  The 
other  profiles  are  approximated  in  the  same  fashion,  with 
the  results  included  in  Table  4.7. 

(1)  Step  One.  As  visibility  decreases  the 
B-coeff icient  becomes  more  negative,  causing  transmittance 
to  decrease  at  a  faster  rate  (more  extinction).  Testing  of 
various  functions  was  conducted  to  generate  the  B— 
coefficient  as  a  function  of  visibility.  It  was  assumed  to 
start  with  that  wind  is  zero  so  that  the  curve  fit  process 
uses  the  B-coeff icients  of  the  first  row  on  Table  4.5. 
Graphically  the  data  points  appear  as  in  Figure  4-4. 

Figure  4—4 


B— 1 Coefficient  Curve  fK  to  Eq.  4.9 


Visibility  In  Kltomotors 
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Care  must  be  taken  in  selecting  a  curve— fit 


method  to  approximate  any  of  the  three  functions  cc<V),  f  (W)  , 
or  6  since  the  B-coef f icient  data  points  are  negative.  This 
eliminates  exponential  or  logarithmic  function  Least-Squares 
techniques  since  the  natural  logarithms  of  the  data  points 
are  necessary,  and  are  undefined. 

To  approximate  « (V) ,  a  cubic  function  was 
selected  of  the  form: 

Row  Ea  (Vis)  Is  A,  +  A,»V  +  As-V*  +  A,-V»  (Eq.  4.9). 

This  eliminated  the  negative  data  point  problem  and  produced 
highly  accurate  results  for  the  few  data  points  on  hand. 

The  values  found  for  the  coefficients  of  the  cubic  are 
contained  in  Table  4.7.  These  values  were  used  to  generate 
the  fitted  curve  in  Figure  4-4,  while  the  original  B- 
coefficients  were  plotted  as  the  single  point  values. 
Remember:  fitting  of  this  function  is  performed  only  on  the 

first  row  of  B— coefficients  in  Table  4.5. 

(2)  Step  Two  is  to  establish  a  function  which 
reflects  the  change  in  the  decay  constant  with  wind:  i.e. 
f <M) .  Graphically  the  data  points  plot  as  shown  in  Figure 
4-5.  At  first  the  function  appears  as  an  exponential  decay 
which  would  be  easily  modelled.  But,  since  the  data  points 
are  negative  this  method  is  not  easily  performed.  This  is 
where  the  concept  of  Table  4.6  is  useful.  Given  that  the  B- 


V,  A 
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coefficient  for  any  visibility  condition  (with  wind  zero) 
can  be  generated  by  completing  step  one,  then  f (W>  becomes 
merely  a  multiplier  to  that  value  generated. 

To  find  out  what  to  multiply  « (V)  by  to  get 
the  decay  constant  for  wind  conditions  which  are  not  zero, 
we  fit  a  function  to  the  first  col umn  of  Table  4.6.  Now, 
all  the  data  points  are  positive,  and  any  function  is  fair 
game.  Graphically,  the  newly  translated  function  appears  as 
in  Figure  4-6.  After  unsuccessful  (inaccurate)  attempts  to 
fit  a  hyperbolic  tangent  (tanh)  function  to  this  curve,  a 
much  more  accurate  one  was  discovered: 

Column  Cf (Wind) 3 :  1  +  ED,  •(  1-Exp  (D2- Wind)  )  3  (Eq.  4.10) 

(3)  Step  Three.  Figure  4-6  shows  a  family  of 
five  curves  which  plot  each  of  the  five  columns  of  Table 
4.6.  The  curve  with  the  least  maximum  represents  the  first 
column  in  Table  4.6.  The  others  in  increasing  fashion 
represent  decreasing  visibility  conditions.  The  correlation 
factor  6  is  necessary  to  increase  the  extinction  when  wind 
and  visibility  are  not  zero  and  maximum  respectively.  In 
this  case  6  acts  as  a  multiplier  also,  driving  the  decay 
constant  more  negative  when  its  value  becomes  greater  than 
one.  Fitting  a  function  to  points  vertically  spaced  on 
Figure  4-6  is  the  key  to  finding  the  correlation  coefficient 
6.  A  function  which  works  extremely  well  in  this  case  is: 


<5  =  C,-Vc* 


(Eq.  4.11) 


The  same  effect  may  be  obtained  by  -fitting 
equation  4.11  to  the  data  points  o-f  one  o-f  the  lower  rows  o-f 
Table  4.6  if.  they  have  been  normalized  a  second  time  in  the 
horizontal  direction  (i.e.  divide  each  coefficient  in  one 
lower  row  of  Table  4.6  by  the  coefficient  in  the  first 
column  of  that  row). 

This  procedure  completes  the  curve  fitting 
process  for  each  atmospheric  profile,  leaving  eight  constant 
coefficients  which,  when  incorporated  into  equation  4.3, 
yield  an  admirable  approximation  of  transmittance  values  at 
range  values  between  one  and  twenty  kilometers  throughout 
the  entire  atmospheric  profile.  This  procedure  was  used  to 
establish  the  constants  for  each  of  the  six  LOWTRAN 
profiles.  These  constants  are  included  here  for  reference. 

TABLE  4.7 


CURVE  FIT  COEFFICIENTS  FOR  ALL  SIX  LOWTRAN  PROFILES 


TROPICAL 

MIDLATITUDE 

MIDLATITUDE 

PROFILE 

SUMMER 

WINTER 

A, 

-0. 139154 

-0. 1249852 

-0.08956 

Aj 

6. 2887E-03 

6. 296102E-03 

6.31731 7E— 03 

A, 

-2. 2354E-04 

—2. 207060E— 04 

-2 . 247950E— 04 

A, 

2. 6091E-06 

2.61 9482E-06 

2. 627775E-06 

c, 

13.74105 

15. 18516 

29.876 

c. 

-0. 716806 

-0. 743405 

-0.923857 

D, 

0.895 

1.10 

2.46 

D, 

-0.30 

-0.27 

-0.295 

or 

0.880 

0.875 

0.9383 

TABLE  4.7  (Continued) 


SUBARCTIC 

SUBARCTIC 

STANDARD 

SUMMER 

WINTER 

ATMOSPHERE 

A, 

-O. 109255 

-0.07970 

-0.091037 

A, 

6.2545E-03 

6. 4225E-03 

5. 6255E-03 

A, 

-2.2228E-04 

—2. 2857E— 04 

—2. 002 IE— 04 

A« 

2.5953E-06 

2.6723E-06 

2. 3407E-06 

c, 

19. 1121 

37.730 

24. 130 

Ca 

-0.80492 

-0.98561 

-0.86706 

D, 

1.440 

3.980 

1.812 

Da 

-0.30 

-0.2825 

-0.34 

0. 

0.9636 

0.9208 

Note  here  that  eight  of  the  nine  constants 
in  each  atmospheric  profile  eventually  constitute  «,  the  'B' 
or  extinction  coefficient.  The  next  step  is  to  implement 
and  explain  the  computer  coding  process.  This  procedure  is 
contained  in  appendix  B.  We  can  now  turn  our  attention  to 
the  SIREM  model  and  examine  the  acquisition  process. 


V.  THE  STOCHASTIC  INFRARED  ENGAGEMENT  MODEL  (SIREM) 


As  discussed  in  chapter  one ,  the  most  inexpensive  means 
of  developing  equipment  characteristics  and  tactics  for 
their  use  is  through  simulation.  The  use  of  computers  to 
analyze  effectiveness  of  systems  or  subsystems  is  becoming 
the  predominant  development  method  for  this  reason.  The 
SIREM  model  was  designed  to  investigate  the  effectiveness  of 
deploying  decoys  which  emit  an  infrared  signature  against 
infrared  seeking  anti-ship  missiles  (ASMs). 

By  examining  through  simulation  such  parameters  as  burn 
time,  deployment  altitude,  deployment  angle,  decoy  drift  and 
descent  characteristics  if  airborne,  and  distance  from  the 
ship,  optimum  design  requirements  and  deployment  tactics  can 
be  developed.  As  with  any  software  model ,  to  maintain  its 
validity  it  must  be  periodically  reviewed  and  upgraded  to 
remove  any  synthetic  qualities  and  replace  them  with 
accurate  representations. 

There  are  ongoing  efforts  to  incorporate  atmospheric 
losses,  discriminative  seeker  capabilities,  and  last  but  not 
least  the  multi-ship  mul ti -mi ssi 1 e  analysis  capability. 
Coincident  with  the  work  being  done  on  the  SIREM  model, 
extensive  investigation  is  being  conducted  regarding  optimum 
deployment  characteri sti cs  for  chaff,  repeater  decoys,  and  a 
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variety  of  other  off board  countermeasures  in  order  to 
integrate  the  use  of  off board  devices  into  a  multi -band 
defensive  system.  To  fully  understand  the  SIREM  model 
requires  extensive  programming  experience  and/or  simulation 
expertise.  Likewise  an  operational  knowledge  of  available 
devices  is  advisable;  however,  a  general  description  of  the 
program  is  in  order. 

A.  THE  MODEL 

The  SIREM  model  calculates  probabilities  associated  with 
distraction  of  an  IR  homing  missile  by  a  pattern  of  decoys 
launched  by  the  attacked  ship.  Rather  then  the  Monte  Carlo 
method  of  simulation,  a  stochastic  approach  was  utilized  in 
order  to  significantly  reduce  computation  time.  The  program 
uses  realistic  input  parameters  pertinent  to  the  scenario 
under  study,  and  evaluates  statistical  information  about  the 
distraction  effectiveness  which  may  be  output  in  graphical 
or  numeric  form. 

The  user  may  easily  modify  any  of  the  input  parameters 
under  the  following  general  headings: 

1.  Ship  characteri st i cs 

2.  Missile  characteristics 

3.  Decoy  characteristics 

4.  Scenario  and  run  characteristics 


Although  the  original  version  transported  -from  Naval 
Research  Laboratory  Mas  designed  to  operate  under  the  NLX 
NAMELIST  -facility  o-f  the  Prime  (TM)  computer  there, 
significant  changes  were  required  to  adapt  the  input/output 
to  non— NAMELIST  type  since  it  was  not  available  at  Naval 
Postgraduate  School;  and  to  modify  the  graphics  subroutines 
and  hardware  to  operate  on  the  Digital  VAX  11/780  (TM)  and 
the  DI— 3000  (TM)  graphics  software  available  there. 

Output  from  the  program  is  available  in  the  form  of 
graphs  or  printed  output.  A  summary  of  the  available  graphs 
is  listed  here: 


1.  Probability  of  distraction  conditioned  on  initial 
mi ssi le-to-ship  distance. 

2.  Probability  of  distraction  conditioned  on  initial 

azimuthal  approach  angle. 

3.  Pattern  of  ship  lock— on  range  as  a  function  of  angle. 

4.  Probability  of  distraction  conditioned  on  initial 
azimuthal  approach  angle  and  a  single  initial  missile 
to  ship  distance. 

5.  Angular  regions  where  probability  of  distraction 
exceeds  a  threshold. 

6.  Locations  where  threshold  crossings  of  plot  5  occur. 

7.  Probability  of  distraction  conditioned  on  initial 
mi ssi le-to-ship  distance  and  sector  of  initial 
azimuthal  approach  angle. 

8.  Ship  orientation,  decoy  deployment,  and  wind 
direction. 

9.  Effects  of  atmospheric  refraction  on  IR  at  various 
incidence  angles. 
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Being  stochastic  by  design,  certain  parameters  may  be 
considered  random  variables  -for  initialization.  The  initial 
position  and  orientation  o-f  the  missile  relative  to  the  ship 
is  a  good  case.  It  is  assumed  -for  startup  that  the  missile 
is  pointed  radially  inward  toward  the  ship,  and  hence  the 
effectiveness  computations  would  yield  a  lower  bound  for 
effectiveness  of  decoy  distraction.  This  is  the  basis  for 
the  'worst  case'  situation  described  in  earlier  chapters. 

The  initial  range  and  azimuth  approach  angle  of  the 
missile  are  discrete  random  variables.  The  initial  range 
may  be  changed  to  appropriate  kilometer  values;  while  the 
initial  range  may  also  be  modified,  it  is  designed  to  run 
multiple  missiles  at  the  ship  every  360/NThct*  degrees. 

The  seeker  characteristics  and  the  random  variables 
assigned  to  emulate  it  can  change  from  missile  to  missile. 

In  this  sense,  randomness  must  be  assigned  to  accomodate 
scanning  directions,  boresight  error,  field  of  view,  and 
probability  of  acquisition  on  multiple  targets  (i.e.  ship 
plus  any  decoys).  The  model  performs  these  adequately, 
using  significant  but  valid  assumptions  regarding  the 
process. 

The  maximum  initial  range  of  the  missile  to  the  target 
is  R.  at  the  time  of  launch  of  the  first  decoy.  This  is 
considered  time  t=0.  It  then  takes  R»/V*  seconds  for  the 
missile  to  reach  the  ship  at  the  origin,  with  V*  being  the 
missile  velocity.  The  sample  time  T4  interval  determines 


the  maximum  number  o-f  time  samples  taken  Nt« .  T,  is  set 


small  enough  to  obtain  an  accurate  representation  of  the 
scenario.  Since  the  missile  is  always  aimed  at  ship  before 
acquisition,  a  lower  bound  on  distraction  results,  and  if  no 
decoy  is  acquired  then  the  probability  of  distraction  is 
zero  for  that  given  initial  range  and  azimuth. 

The  choice  of  time  interval  Tt  is  dependent  on  the 
missile  velocity,  target  velocity,  and  the  scanner  sweep 
time  T».  Since  the  area  of  the  scanner  field  of  view  on  the 
sea  surface  is  large  relative  to  the  distance  either  the 
missile  or  ship  travels  during  the  sweep  time,  Ta  may  equal 
T, .  As  different  conditions  arise  a  necessity  for  time 
scaling  may  result.  Developments  have  begun  to  incorporate 
different  scanner  parameters. 

The  projection  of  field  of  view  is  incorporated  to  allow 
both  curved  earth  and  flat  earth  calculation.  The  ship  and 
decoys  are  treated  as  point  targets  for  infrared  radiator 
purposes.  The  "hot  spot"  of  the  deployed  decoy  is  placed  at 
a  known  distance  above  the  sea  surface.  Calculations  are 
made  to  determine  if  the  hot  spot  is  within  the  field  of 
view  of  the  scanner.  The  next  stage  of  development  requires 
incorporating  radiance  and  transmittance  calculation,  and 
detector  sensitivity  parameters  to  determine  if  there  is 
sufficient  radiance  from  the  ship  or  decoy  to  be  detected  by 
the  seeker.  Much  of  this  development  was  covered  in  earlier 
chapters. 
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The  model  can  also  calculate  any  shadowing  o-f  decoys  by 
the  ship.  This  means  that  if  the  aspect  o-f  the  ship  comes 
between  the  missile  and  decoy,  effectiveness  of  that  decoy 
will  inherently  drop  to  zero.  This  can  be  done  using  either 
the  flat  or  4/3  earth  model.  In  addition  to  shadowing  by 
the  ship,  if  excessive  range  exists,  shadowing  by  the  earth 
may  occur. 

Once  acquisition  of  a  decoy  occurs,  the  missile  flies 
according  to  the  aerodynamic  parameters  of  the  control 
surfaces  and  systems  for  whatever  missile  is  being  modelled. 
If  a  decoy  is  not  in  the  field  of  view  and  the  ship  is 
acquired,  the  missile  is  assumed  to  hit  the  ship.  If  locked 
onto  a  decoy,  the  missile  flies  at  the  decoy  until  it  passes 
it.  At  that  point  break-lock  occurs.  This  condition  may 
also  occur  if  the  control  system  forces  the  missile  to 
maneuver  so  that  the  decoy  is  no  longer  in  the  field  of 
view.  And  lastly,  if  the  decoy  burns  out  break— lock  occurs. 

When  break-lock  occurs  the  scanner  must  re-enter  the 
acquisition  mode.  Since  the  missile  usually  makes  drastic 
maneuvers  near  break— lock,  the  re-acquisition  process  is 
somewhat  random,  but  if  the  ship  is  in  the  field  of  view  it 
may  be  assumed  to  hit  the  ship.  Missile  flight  continues  at 
break-lock  until  all  targets  are  passed. 

The  flexibility  of  the  program  is  excellent.  Virtually 
any  parameter  is  changeable  either  at  input  or  at  source 
code  level.  As  many  as  five  decoys  may  be  deployed  by  the 


ship  (more  if  certain  array  dimensions  are  changed)  in  any 
combination  o-f  deployment  angles,  times,  altitudes,  and 
active  li-fetime.  The  ship  may  change  size  and  speed.  The 
wind  is  a  -factor  on  decoy  drift  and  missile  flight.  And  the 
missile  may  change  speed,  altitude  of  approach,  and  control 
system  or  seeker  parameters.  The  variations  are  endless, 
and  analysis  of  the  many  combinations  pertinent  to  real- 
world  scenarios  may  be  undertaken  to  make  maximum  effective 
use  of  off board  countermeasures.  Additionally,  by  further 
extending  not  only  SIREM  but  other  similar  models  as  well, 
evaluation  of  descrimination  techniques  and  optimization  of 
decoy  signatures  become  possible  avenues  of  further 
development.  With  this  in  mind  we  turn  our  attention  to  the 
acquisition  process  and  observed  phenomena. 


V I •  IR  SEEKER  SYSTEM  AND  COUNTERMEASURES 

Since  the  principal  interest  here  lies  in  passive 
seekers,  it  is  convenient  to  separate  these  into  three  broad 
categories:  (1)  Scanning  systems  such  as  search,  track,  and 
imagers,  (2)  Staring  systems  which  use  non-mechanical  means 
to  obtain  spatial  scene  viewing,  and  (3)  Measurement  systems 
such  as  radiometers,  spectrometers  and  interferometers.  Of 
interest  here  are  the  first  two  categories.  Scanning 
detection  schemes  are  most  common  since  staring  systems 
usually  require  sophisticated  mul t i -di mensi onal  array 
detectors  and  massive  numbers  of  amplifiers  to  support  the 
individual  elements. 

A.  SCANNING  AND  STARING  SYSTEMS 

Scanning  systems  sample  the  radiant  intensity 
distribution  in  their  field  of  view  using  various 
hori zontal /vert i cal  techniques.  The  output  from  the 
detector  in  this  case  will  be  a  linear  analog  of  the 
distribution  or  a  simple  indication  of  the  presence  of  a 
target  and  its  location.  The  object  plane  is  usually 
assumed  at  infinity  and  the  detector  responds  to  radiant 
levels  within  the  instantaneous  field  of  view  of  its 
element (s).  This  gradient  in  radiance  is  translated  to 
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electrical  signals  which  are  then  filtered,  amplified,  and 
referenced  to  a  baseline  (missile  axis)  for  inputs  to  the 
various  azimuth  and  elevation  control  systems.  For  imaging 
systems  the  signal  is  translated  from  the  object  to  its 
corresponding  image  plane  by  painting  out  the  radiance 
distribution  using  optical  deflection  onto  a  film,  or  beam 
deflection  to  a  cathode  ray  tube.  Using  internal  processing 
techniques  in  missile  seekers,  the  translation  is  used  as 
input  to  the  missile  control  surface  system  to  center  the 
highest  radiant  intensity  in  the  center  of  the  total  field 
of  view.  The  use  of  either  spiral  or  raster  techniques  to 
accomplish  slow  mechanical  scanning  is  most  common.  In  some 
cases  the  relatively  high  speed  motion  at  altitude  provides 
the  vertical  trace  capability  for  the  system,  while  a 
spinning  mirror  provides  the  horizontal  trace.  Regardless 
of  the  scanning  technique  employed,  major  concern  has 
recently  been  given  to  the  incorporation  of  discrimination 
processing  into  newer  seeker  systems.  If  (or,  more 
appropriately,  when)  this  occurs,  the  effectiveness  of  using 
offboard  decoys  could  be  significantly  reduced.  Therefore, 
an  examination  of  possible  discrimination  methods  is  in 
order . 

In  dealing  with  imaging  systems  there  are  two  basic 
approaches  to  describe  the  sensitivity  of  the  detector.  The 
first  approach  is  based  on  the  detector  figure  of  merit  (D~> 
described  earlier.  It  regards  the  detector  as  a  noise 


source  of  -fixed  level  which  is  determined  from  the 


measurements  of  fixed  conditions  similar  to  those 
encountered  by  the  detector  when  installed  in  its  normal 
operating  configuration  inside  the  seeker. 

The  other  approach  which  is  suitable  only  for  BLIP 
(Background  limited  mode)  detectors  regards  a  random  arrival 
rate  of  photons  from  the  background  as  the  ultimate  source 
of  detector  noise.  In  such  a  case  the  detector  noise  power 
density  is  described  using  the  laws  of  photon  rate  and 
detector  quantum  efficiency.  When  the  number  of  incoming 
(target)  photons  is  small  in  comparison  to  the  number  of 
background  photons  the  approaches  are  equivalent,  since  BLIP 
D*  is  determined  by  background  photon  rate.  If  the  number 
of  target  photons  is  larger  than  the  average  number  of 
background  photons  though,  the  noise  due  to  the  random 
arrival  time  of  target  photons  dominates  and  D*  then  is 
meaningless.  This  case  is  the  unusual  one. 

Image  forming  scanners  are  not  usually  designed  to 
respond  to  an  absolute  radiance  level,  i.e.  not  for  a  dc 
response.  They  are  designed  to  follow  changes  in  radiance 
as  the  field  of  view  is  scanned.  This  includes  gradual 
changes  in  the  radiance  level  corresponding  to  low 
frequencies.  DC  restoration  can  be  used  to  control  the 
intensity  levels  at  the  output.  Without  this  form  of 
automatic  gain  control  (AGO  periods  of  signal  blackout  or 
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saturation  would  occur  due  to  high  l/-f  noise  and  amplifier 
drift  at  low  frequencies. 

The  signal  of  interest  from  the  scene  is  the  variation 
in  irradiance  H  caused  by  the  variation  in  temperature  or 
emissivity  of  objects  in  the  scene  and  may  be  expressed  as: 

a  a> 

AH  =  [l  /at f^M  d\  +  f A€i.  WvdXl  in  Matts  (Eq.  6.1) 

\  J  JT  J  €„  /  cm* 

o  o 

where  V  -  instantaneous  field  of  view  in  steradians 
T  =  temperature  in  degrees  Kelvin 
=  emissivity  (wavelength  dependent) 

Ux  =  spectral  radiant  emittance  in  M/cm* 

If  most  noise  in  the  system  is  caused  by  the  detector, 
it  may  be  expressed  in  terms  of  "noise  equivalent  flux 
density"  (NEFD) .  This  is  the  minimum  signal  irradiance 
capable  of  producing  a  peak  signal  to  rms  noise  ratio  (Sp/Nn) 
of  one. 

The  ultimate  goal  is  to  find  the  actual  signal-to-noise 
ratio  (S/N)  generated  by  the  scanner.  This  is  given  as: 

S  =  Variation  of  irradiance  in  scene  (AH)  in  Matts  (Eq.  6.2) 
N  NEFD 

The  noise  equivalent  power  (NEP)  in  watts  of  the  system 
can  be  obtained  from  the  figure  of  merit  of  the  detector 
(D~)  and  is  expressed  as: 
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(Eq.  6.3) 
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where  A*  =  area  of  the  detector  in  cm2 
Af  =  system  noise  bandwidth 

D**  =  detector  figure  of  merit  in  cm-Hzl/2/watt 

The  signal  power  on  the  detector  is  equal  to  the  optical 
collector  area  times  the  signal  flux  density.  Therefore  the 
system  noise  equivalent  flux  density  is: 

NEFD  =  4  Vfo-  ■  Af  in  Watts/cm2  (Eq.  6.4) 

*D.2D-€. 

where  D„  =  diameter  of  the  aperture  of  the  optical  system 
€.  =  optical  efficiency 

The  difference  in  temperature  and  emissivities  (AT  and 
A€)  of  objects  in  the  scene  contribute  to  the  differential 
radiant  emittance  (AH),  but  due  to  the  small  contribution  of 
A€  for  non-BLIP  detector  operation,  its  term  in  equation  6.1 
is  commonly  neglected.  By  substituting  W*  for  the  first 
integral  term  the  abbreviated  value  for  radiant  emittance 
becomes  then: 

AH  =  V*ATWf  in  Watts/cm2  (Eq.  6.5) 

s’ 

The  instantaneous  field  of  view  v  is  defined  by  the 
detector  area  Ad  and  the  focal  length  fl  as: 

V  =  A*  in  steradians  (Eq.  6.6) 

(f  l)2 
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IT: 


By  substituting  equation  6.6  into  equation  6.4  the  noise- 


equivalent  flux  density  is  reestablished  as: 


NEFD  =  4F  V»-Af 
*D.D~€. 


in  Watts/cm2 


(Eq.  6.7) 


where  F  is  the  f-number  of  the  optics. 

The  signal-to-noise  ratio  generated  by  the  scanning 
system  is  then  (from  equation  6.2): 


S  =  &TH*D.D<>€.yi» 
N  4F  Vfff 


(Eq.  6.8) 


Some  careful  considerations  are  in  order  here  since  the 
AH  actually  seen  by  the  sensor  is  modified  by  the 
transmittance  of  the  optical  components  and  the  atmosphere 
as  discussed  earlier.  If  the  bandwidth  considered  is 
sufficiently  narrow,  transmittance  losses  may  be  compensated 
for  by  using  the  average  transmittance  values  for  each  while 
still  producing  accurate  results.  The  narrow  bandwidth  was 
the  premise  for  the  development  of  transmittance  modelling 
in  Chapter  IV. 

Sensitivity  of  an  image  forming  scanner  is  expressed  by 
its  noise  equivalent  temperature  (NET).  This  parameter 
refers  to  the  temperature  difference  AT  within  a  scene 
element  (with  a  given  emissivity  and  average  temperature  of 
e.g.  300°K)  which  can  produce  a  change  in  the  electrical 
signal  level  equal  to  the  rms  system  noise.  The  noise  in  a 
well  designed  system  is  mostly  detector  noise.  By 
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incorporating  the  atmospheric  transmittance  into  equation 
6.8  and  solving  -For  AT,  the  NET  is  obtained  as: 

NET  =  4F  l/Af~  in  «*K  (Eq.  6.9) 

*D.D-€.T»W, 

This  an  expression  of  the  sensitivity  of  the  detector  system 
as  a  whole.  It  makes  the  assumption  that  the  scene  area 
over  which  the  temperature  difference  AT  occurs  is  much 
larger  than  the  field  of  view. 

From  the  design  aspect  a  major  emphasis  must  be  placed 
on  the  value  of  D~  for  a  given  detector  since  its  value  may 
vary  as  much  as  four  orders  of  magnitude  among  common 
detector  materials  CRef  2:  Table  7.11.  In  addition,  the  D* 
of  a  particular  detector  will  vary  considerably  from  the 
published  data  owing  to  such  factors  as  manufacturing 
process,  electrical  bias,  cooling,  loading,  aging  and 
aerodynamic  heating  of  the  optical  window. 

The  electrical  bandwidth  <Af)  of  the  system  is 
determined  by  the  scan  rate.  It  is  proportional  to  the 
reciprocal  of  the  dwell  time  of  an  image  point  on  the 
scanning  detector  surface.  The  amplifiers  connected  to  the 
detector  can  be  optimized  to  obtain  a  maximum  signal-to- 
noise  ratio,  especially  if  the  characteristics  of  the  signal 
and  noise  are  known.  In  this  fashion  a  number  of 
target/ob ject  discrimination  techniques  may  be  devised. 
Optimum  filter  design  is  extremely  difficult  to  perform 
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using  only  a  -few  components,  though  a  simple  band-pass 
filter  with  a  flat  response  and  a  low  and  high  frequency 
cutoff  is  often  adequate.  The  very  high  signal  to  noise 
ratios  realized  by  using  matched  filter  (correlation) 
techniques  for  processing  are  virtually  impossible  due  to 
the  passive  nature  of  the  systems,  but  by  tweaking  the  high 
and  low  frequency  roll-off  a  much  improved  background 
control  and  target  definition  may  be  accomplished  with  only 
a  slight  reduction  of  the  S/N  ratio. 

If  a  total  field  of  view  of  fl  steradians  is  scanned  at  a 
constant  rate  without  redundancy  and  with  an  instantaneous 
field  of  view  W  steradians,  then  the  scan  is  completed  in  a 
frame  time  of  ta.  The  dwell  time  t*  of  a  point  on  the 
detector  is  them  As  the  total  field  of  view  is 

increased  (either  in  azimuth,  elevation  or  both)  linearity 
in  scan  rate  is  difficult  to  achieve  CRef  4:pp  7351.  The 
off-center  loss  in  dwell  time  must  be  compensated  for  by  a 
scan  efficiency  factor  €«.  As  the  scan  efficiency 
decreases,  an  increased  electrical  bandwidth  is  required. 

In  this  case  the  optimum  bandwidth  required  is 
approximately: 

Af„  =  3Q€.  in  Hz  (Eq.  6.10) 

4*t. 

Actually  this  refers  to  the  upper  cutoff  frequency  since  it 
is  not  practical  to  design  for  dc  response.  The  lower 


cutoff  frequency  can  usually  be  designed  to  handle  down  to 
few  Hertz  which  allows  reproduction  of  at  least  most  of  the 
low  frequency  features  in  the  scene.  If  a  linear  (vertical 
usually)  array  of  detectors  is  used,  with  N  elements  in  the 
array,  then  Af„  can  be  reduced  proportionately . 

Equivalent  noise  bandwidth  Af  will  usually  be  larger 
than  required  bandwidth  Af«  because  of  the  inherent  1/f 
noise,  amplifier  noise,  and  R-C  filter  roll-off.  In  terms 
of  scan  rate,  the  noise  equivalent  temperature  of  a  scan 
system  is: 


2  3€.  in  -K  (Eq.  6.11) 
€.T. 


By  breaking  equation  6.11  into  the  above  groups  the 
corresponding  components  may  be  individually  analyzed.  The 
first  group  represents  system  performance  requirements.  Th 
second  group  represents  mostly  the  physical  components 
available  at  design  inception,  and  the  third  group 
represents  efficiency  parameters  of  materials  used  and  the 
atmosphere. 

The  principles  of  detectors  and  line  scanning  systems 
are  covered  extensively  in  References  C21  and  C33.  The 
basics  presented  thus  far  will  serve  as  a  conceptual 
approach  to  modelling  seeker  systems,  and  will  serve  as  a 
baseline  for  devising  discrimination  techniques  in  the 
modelling  process.  Based  on  discussion  thus  far  no 


committment  has  been  made  with  regards  to  any  particular 
system,  but  either  conceived  or  existing  seekers  which  use 
the  line  scanner  techniques  operate  under  the  previously 
explained  principles. 

1 .  Filtering  for  Discrimination 

There  are  three  -filtering  methods  which  may  be  used 
in  seeker  system  designs.  These  are:  1)  Spectral,  2) 
Spatial,  and  3)  Temporal.  Each  method  has  its  physical 
limitations  -for  implementation  in  discrimination  as  will  be 
discussed  shortly.  In  contrast,  the  incorporation  o-f 
discrimination  into  IR  seekers  leads  to  the  necessity  to 
analyze  present  countermeasures  for  possible  enhancement  to 
subjugate  any  of  these  discrimination  techniques. 

a.  Spectral  Filtering 

Since  background  radiant  intensity  is  a  primary 
consideration  to  seeker  design,  methods  to  reduce  it  are  a 
necessity.  Background  radiant  intensity  typically  has  peaks 
(maxima)  at  approx imately  0.5  and  10  micrometers  caused  by 
solar  reflection  and  thermal  self-emission  respectively. 
Around  3  microns  a  region  of  minimum  background  radiance 
exists.  Coincidentally,  between  the  two  maxima,  most  hot 
targets  produce  their  peak  radiant  intensity.  The  simplest 
method  of  detection  is  to  examine  typical  target  spectral 


characteri  sties,  select  a  detector  most  appropriate  -for 
operation  at  the  target's  peak  spectral  points,  and  then  use 
filtering  to  prevent  the  background  radiant  energy  from 
entering  the  system.  If  the  target  of  interest  exhibits 
significant  wavelength-dependent  emissions,  one  method  of 
discrimination  could  be  the  use  of  notch  (2-color  spectral 
bandpass)  filters  to  select  specific  wavebands.  If  the 
radiant  energy  in  the  passbands  produced  by  the  decoy  is 
then  less  than  that  produced  by  the  target  (W/sr) ,  the  decoy 
may  be  rendered  ineffective  in  the  seduction  role.  Thus 
much  of  the  design  of  a  decoy  revolves  around  matching  its 
spectral  radiance  to  that  of  the  ship  it  is  being  designed 
to  protect. 

Most  IR  decoys  are  omni-directional  by  nature, 
and  exhibit  the  same  spectral  signature  in  all  directions 
(disregarding  wind  or  sea-surface  effects  if  floating). 
Target  signatures  are  aspect  dependent  (intensity)  as  well 
as  spectrally  dependent.  As  the  aspect  changes,  so  does  the 
intensity  of  certain  spectral  components.  The  inherent 
variability  of  target  signatures  in  this  respect  makes  it 
exceedingly  difficult  to  establish  discrimination  techniques 
since  look-up  tables  and  matched  filter  methods  cannot  be 


b.  Temporal  Filtering 


Spatial  frequencies  produce  time  varying  signals 
with  a  frequency  content  equal  to  the  spatial  frequency 
times  the  angular  scan  rate.  In  processing  then,  low 
frequency  rejection  and  high  frequency  peaking  can  be  used 
to  eliminate  background  clutter.  Only  if  the  target  and 
background  signals  can  be  adequately  described  can  a  spatial 
filter  maximize  the  ratio  of  target  signal  to  background 
signal.  If  mul ti -el ement  detectors  are  used,  simpler  filter 
techniques  are  possible.  The  instantaneous  field  would  be? 
made  as  small  as  practical ,  and  the  upper  frequency  cutoff 
of  the  electrical  filter  is  selected  to  correspond  to  the 
dwell  time  of  a  point  on  the  detector.  The  low  frequency 
cutoff  is  usually  set  at  1/3  to  1/4  of  the  upper  cutoff 
frequency  to  block  the  background  radiance  without  allowing 
saturation  at  the  other  end  of  the  bandpass.  This  filtering 
technique  is  most  often  used  to  decrease  the  probability  of 
false  alarms  through  scan-to— scan  comparisons  and  has  little 
merit  in  discerning  targets  from  decoys. 

c.  Spatial  Filtering 

A  more  promising  approach  to  discrimination  by 
the  seeker  may  be  to  process  spatial  character i sties  of 
scene.  When  the  range  from  the  seeker  to  the  target  is 
large,  the  objects  in  the  scene  will  appear  as  point 
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sources.  As  this  range  decreases,  however,  large  objects 
form  a  distribution  within  the  scene  (both  horizontal  and 
vertical)  while  small  objects  retain  their  point  source 
characteristics.  By  keeping  the  instantaneous  field  of  view 
small,  integral  (in  spatial  sense)  processing  could  be  used 
to  discern  large  distributed  objects  from  small  ones.  It 
then  becomes  the  task  of  the  decoy  designer  to  not  only  make 
the  decoy  more  intense  spectrally  than  the  ship,  but  to  make 
the  signature  more  widely  distributed  to  simulate  large 
targets.  If  the  decoy  is  suspended  in  the  air,  for  example 
by  parachute  (and  somewhat  less  so  if  the  decoy  is 
waterborne),  the  wind  will  tend  to  distribute  the  heat  over 
a  wider  area;  but  variable  wind  conditions  and  short  hang 
times  of  parachute-deployed  decoys  make  this  method  of 
deployment  less  attractive  than  water — borne  placement.  Some 
of  the  effects  of  wind  for  signature  expansion  are  still 
present  with  water-borne  decoys.  In  addition  to  the  torch 
emission  above  the  sea  surface  there  exists  a  reflective 
component  off  the  sea  surface  referred  to  as  the  "Halo 
effect"  CRef.  43.  This  effect  was  originally  noticed  by 
personnel  from  Naval  Research  Laboratory  taking  infrared 
measurements  of  the  USS  Ti conderoga ' s  signature  in  both  the 
3-5  and  8—14  micron  bands.  The  halo  effect  may  readily  be 
seen  in  Figure  6-1.  Upon  investigation  it  was  established 
that  the  sea  foam  caused  by  a  white-capped  sea  surface  has 
very  high  reflectance  at  low  grazing  angles  in  both  infrared 
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bands.  The  white-caps  may  be  caused  by  wind,  ship  movement, 
and  other  surface  turbulence.  When  viewed  from  certain 
angles  as  in  Figure  6-2,  energy  from  the  source  which  was 
not  originally  emitted  in  the  direction  of  the  seeker 
(observer)  is  reflected  from  the  sea  surface  toward  the 
seeker,  adding  to  the  original  signature  radiant  density. 

From  the  decoy  design  aspect  this  effect  may  be 
capitalized  upon  by  i ncorporati ng  a  means  of  disturbing 
(foaming)  the  sea  surface  over  a  wide  area  around  a  water — 
borne  decoy  beneath  its  burning  source.  By  floating  the 
decoy  at  the  proper  height  above  the  foamed  sea  surface  a 
very  large  reflective  signature  may  be  achieved  in  addition 
to  the  intense  point  source  of  the  torch  itself. 

Foaming  may  be  accomplished  in  many  ways,  but 
primary  interest  must  be  in  distributing  over  as  large  an 
area  as  possible.  Since  it  is  difficult  to  foam  the  sea 
surface  over  a  wide  area  using  a  simple  squi rrel -cage  fan  or 
compressed  air  (too  many  cylinders)  another  prospect  must  be 
taken  which  considers  weight  and  space  restrictions  of  the 
decoy.  Foaming  of  sea  water  is  purely  transient  and  very 
short  in  duration,  so  the  above  methods  would  probably  be 
discarded  upon  testing.  The  alternate  means  suggested  here 
is  to  use  a  small  pneumatically  driven  suction  pump  with  a 
mixing  chamber  to  mix  the  incoming  sea  water  with,  for 
example,  light-water  or  aqueous  film-forming  foam  (AFFF) . 

The  mixture  may  then  be  sprayed  from  the  decoy  by  fan  or 


nozzles  over  a  wide  area  of  the  sea  surface  to  create  a 
highly  reflective  layer  (using  sea  surface  turbulence  for 
wide  angular  reflection)  and  correspondingly  large 
signature.  The  components  in  this  case  could  be  very 
inexpensive.  Additionally,  this  technique  would  work  in  all 
infrared  bands  and,  because  of  the  difficulty  in  devising 
passive-only  discrimination  methods,  may  be  the  optimum 
off board  countermeasure  to  IR  seekers.  The  methods 
discussed  here  have  not  been  attempted,  but  these  and  many 
other  approaches  will  need  to  be  addressed  in  order  to  meet 
the  IR  threat  with  inexpensive  off board  countermeasures. 


-  TABL8  A.l.a  - 

TRANSMITTANCE  IN  3.5  ON  BAND  -  TROPICAL  PROFILE 
***  LOWTRAN  DATA  *** 

Range  !  Visibility  in  Kiloaetera 


in  Ka  ! 

40 

30 

20 

15 

10 

5 

Wind=0 

0.20 

93.17* 

93.14* 

93.07* 

93.01* 

92.88* 

92.50* 

1.00 

82.05* 

81.91* 

81.63* 

81.35* 

80.80* 

79. 10* 

2.00 

73.19* 

72.95* 

72.46* 

71.97* 

70.99* 

68.16* 

4.00 

60.33* 

59.92* 

59. 11* 

58.30* 

56.73* 

52.28* 

8.00 

43.64* 

43.06* 

41.90* 

40.77* 

38.59* 

32.77* 

12.00 

32.70* 

32.09* 

30.81* 

29.58* 

27.27* 

21.40* 

18.00 

22.12* 

21.46* 

20.18* 

18.98* 

16.77* 

11.61* 

21.00 

18.45* 

17.86* 

16.64* 

15.50* 

13.45* 

8.80* 

Wind=5 

0.20 

92.16* 

91.75* 

90.93* 

90.11* 

88.49X 

83.82* 

1.00 

77.70* 

75 . 99* 

72.64* 

69.42* 

63 . 40* 

48.37* 

2.00 

65.70* 

62.79* 

57.39* 

52 . 40* 

43.71* 

25.45* 

4.00 

48.50* 

44.39* 

37.07* 

30.91* 

21.50* 

7.28* 

8.00 

28.30* 

23.63* 

16.48* 

11.46* 

5.55* 

0.64* 

12.00 

17. 10* 

13.05* 

7.60* 

4.40* 

1.50* 

0.06* 

18.00 

8.30* 

5.56* 

2.47* 

1.09* 

0.21* 

0.03* 

21.00 

5.90* 

3.70* 

1.44* 

0.55* 

0.08X 

.00* 

Wind=10 

0.20 

91.76* 

91.47* 

90.49* 

89 . 52* 

87.61* 

82.14* 

1.00 

76.84* 

74.83* 

70.92* 

67 . 19* 

60.31* 

43.70* 

2.00 

64.21* 

60 . 88* 

54.70* 

49.10* 

39.56* 

20.78* 

4.00 

46.41* 

41.73* 

33.68* 

27.13* 

17.61* 

48.60X 

8.00 

25.83* 

20.88* 

13.60* 

8.83* 

3.72* 

0.28* 

12.00 

14.90* 

10.84* 

5.70* 

2.98* 

0.82X 

0.02* 

18.00 

6.80* 

4.21* 

1.61* 

0.61X 

0.09* 

0.01* 

21.00 

4.67* 

2.67* 

0.87* 

0.28* 

0.03* 

.00* 

Wind=15 

0.20 

91.90* 

91.39* 

90.37* 

89.36* 

87.37X 

81.69X 

1.00 

76.61* 

74. 5 IX 

70.46* 

66.59* 

59.49* 

42.51* 

2.00 

63.81* 

60.37X 

53.98* 

48.22* 

38.49* 

19.65* 

4.00 

45.85* 

41.03* 

32.80* 

26. 17* 

16.68* 

4.35* 

8.00 

25.21* 

20.19* 

12.91* 

8.22* 

3.33* 

0.23* 

12.00 

14.38* 

10.30* 

5.27* 

2.68* 

0.69* 

0.00* 

18.00 

6.43* 

3.90* 

1.43* 

0.52* 

0.07* 

0.00* 

21.00 

4.38* 

2.45* 

0.76* 

0.23* 

0.02* 

0.00* 

Wind=2t» 

0.20 

91.87* 

91.35* 

90.32* 

89.28* 

87.26* 

81.47* 

1.00 

76.50* 

74.36* 

70.24* 

66.31* 

59. 11X 

41.96* 

2.00 

63.62* 

60.12* 

53.65* 

47.82* 

38.00* 

19. 15* 

4.00 

45.59* 

40.71* 

32.40* 

25.74* 

16.25* 

4.13* 

8.00 

24.92* 

19.87* 

12.59* 

7.95* 

3.17* 

0.20* 

12.00 

14.08* 

10.03X 

5.08* 

2.55* 

0.64* 

0.01* 

18.00 

6.27* 

3.77* 

1.35* 

0.48* 

0  06* 

0.01* 

21.00 

4.23* 

2.33* 

0.71* 

0.21* 

0.02* 

.00* 
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-  TABLB  A.l.b  - 

TRANSMITTANCE  IN  3.5  UM  BAND  -  TROPICAL  PROFILE 

***  EMPIRICAL  *** 


Range  ! 
in  Kn  !  40 


=  0  0.2 

=  0  1. 

=  0  2. 

=  0  4. 

=  0  8. 

=  0  12. 

=  0  18. 

=  0  21. 


Visibility  in  Kiloaeters 
30  20  15  10 


Wind=  5 

0 

Wind=  5 

1 

Wind=  5 

2 

Wind=  5 

4 

Wind=  5 

8 

Wind=  5 

12 

Wind=  5 

18 

Wind=  5 

21 

Wind=10 

0 

Wind= 10 

1 

Wind=10 

2 

Wind=10 

4 

Wind=10 

8 

tfind=10 

12 

Wind=10 

18 

Wind= 10 

21 

Wind= 15 

0 

Wind= 15 

1 

Nind=15 

2 

Nind=15 

4 

Hind=15 

8 

Nind=15 

12 

Wind=15 

18 

Wind=15 

21 

Wind=20 

0 

W ind=20 

1 

Wind=20 

2 

Wind=20 

4 

Wind=20 

8 

Wind=20 

12 

«ind=20 

18 

Wind=20 

21 

86.63* 

86.58* 

86.57* 

86.49* 

86.33* 

86.03* 

81.37* 

81.13* 

81.08* 

80.72* 

79.94* 

78.60* 

75.25* 

74.80* 

74.70* 

74.05* 

72.63* 

70.20* 

64.34* 

63.59* 

63.41* 

62.31* 

59. 94* 

56.01* 

47.04* 

45.95* 

45.69* 

44.12* 

40.82* 

35 . 64* 

34.40* 

33.20* 

32.93* 

31.24* 

27.80* 

22.68* 

21.50* 

20.39* 

20.14* 

18.61* 

15.63* 

11.52* 

17.00* 

15.98* 

15.75* 

14.36* 

11.72* 

8.21X 

85.75* 

85.14* 

84.16* 

83.06* 

80.76* 

74.54* 

77.30* 

74.60* 

70.42* 

65 . 94* 

57.28* 

38.36* 

67.91* 

63.23* 

56.35* 

49.41* 

37.29* 

16.72* 

52.40* 

45.44* 

36.08* 

27.74* 

15.80* 

3. 18* 

31.21* 

23.46* 

14.79* 

8.74* 

2.84* 

0.11* 

18.58* 

12.11* 

6.07* 

2.76X 

0.51* 

0.00* 

8.54* 

4.49* 

1.59* 

0.49* 

0.04X 

0.00* 

5.79* 

2.74* 

0.82X 

0.21X 

0.01* 

0.00* 

85.55* 

84.88* 

83.82* 

82.63* 

80.13* 

73.41* 

76.39* 

73.48* 

69.00* 

64.22* 

55.08* 

35.56* 

66.32* 

61.35* 

54.10* 

46.87* 

34.47X 

14.37* 

49.98* 

42.77* 

33.26* 

24 . 96* 

13.50* 

2.35* 

28.38* 

20.79* 

12.57* 

7.08* 

2.07* 

0.06* 

16.12* 

10.10* 

4.75* 

2.01* 

0.32X 

0.00* 

6.90X 

3.42* 

1.10* 

0.30* 

0.02* 

0.00* 

4. SIX 

1.99X 

0.53* 

0.12* 

0.00* 

0.00* 

85.50* 

84.82* 

83.74* 

82.53* 

79.99* 

73.16* 

76. 19* 

73.23* 

68.68* 

63.84* 

54.60* 

34.96* 

65 . 97* 

60.94* 

53.61* 

46.32* 

33.87* 

13.89* 

49.45* 

42.19* 

32.66X 

24.38* 

13.04* 

2. 19* 

27.79* 

20.23* 

12.12* 

6.75* 

1.93* 

0.05* 

15.61* 

9.70* 

4.50* 

1.87* 

0.29* 

0.00* 

6.58* 

3.22* 

1.02* 

0.27* 

0.02* 

0.00* 

4.27* 

1.86* 

0.48* 

0. 10* 

0.00* 

0.00* 

85.49* 

84.81* 

83.73* 

82.51* 

79.96* 

73.11* 

76. 15* 

73.17* 

68.61* 

63.76X 

54.49* 

34.83X 

65.89* 

60.84* 

53.50* 

46.19* 

33.74* 

13.78* 

49.33* 

42.07* 

32.52* 

24.25* 

12.94* 

2.16* 

27.66X 

20. 11* 

12.02* 

6.68* 

1.90X 

0.05* 

15.50* 

9.61* 

4.44* 

1.84* 

0.28* 

0.00* 

6.51* 

3.18* 

1.00* 

0.27* 

0.02* 

0.00* 

4.22* 

1.83* 

0.47* 

0.10* 

0.00* 

0.00* 
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-  TABLE  A. 2. a  - 

TRANSMITTANCE  IN  3.5  UM  BAND  -  MIDLATITUDE  SUMMER  PROFILB 
***  LONTRAN  DATA  *** 

Range  !  Visibility  in  Kiloneters 


in  Ka  1 

40 

30 

Wind=0 

0.20  : 

94. 16% 

94.12* 

1.00  : 

84.47* 

84.33* 

2.00  : 

76.79* 

76.53* 

4.00  : 

65 . 59* 

65.16* 

8.oo  : 

50.26* 

49.59* 

12.00  : 

39.82* 

39.03* 

I8.00  : 

28.86* 

28.01* 

21.00  : 

24.90* 

24.04* 

Wind=5 

0.20  : 

93.14* 

92.72* 

l.oo  : 

80.01* 

78.24* 

2.00  : 

68.89* 

65.88* 

4.00  : 

52.78* 

48.28* 

8.00  : 

32.55* 

27.23* 

12.00  : 

20.75X 

15.88* 

18.00  : 

10.86* 

7.27* 

21.00  : 

7.96* 

4.98* 

Wind=10 

0.20  : 

92.93* 

92.44* 

l.oo  : 

79.11* 

77.04* 

2.00  : 

67 . 35* 

63.86* 

4.00  : 

50.46* 

45.37X 

8.oo  : 

29.74* 

24.05* 

12.00  : 

18. 13* 

13.18* 

18.00  : 

8.87* 

5.49* 

21.00  : 

6.28* 

3.60* 

Wind=15 

0.20  : 

92.87* 

92.36* 

l.oo  : 

78.87* 

76.71* 

2.00  : 

66.94* 

63.32* 

4.00  : 

49.84X 

44.60* 

8.oo  : 

29.02* 

23.24* 

12.00  : 

17.47* 

12.52* 

18.00  : 

8.39* 

5.09* 

21.00  : 

5.89* 

3.29* 

Wind-20 

0.20 

92.84* 

92.32* 

1.00 

78.75* 

76.56* 

2.00 

66.74* 

63.07* 

4.00 

49.55* 

44 . 25* 

8.00 

28.68* 

22.87* 

12.00 

17.17* 

12.22* 

18.00 

8.17* 

4.91* 

21.00 

5.71* 

3.15* 

20 

15 

10 

5 

94.06* 

94 . 00* 

93.87* 

93 . 49* 

84.05* 

83.76* 

83.20* 

81.53* 

76.02X 

75.51* 

74.49* 

71.53* 

64.28* 

63.42* 

61.72* 

56.92* 

48.27* 

46.98* 

44.50* 

37.84* 

37.48* 

35 . 99* 

33.17* 

26.02* 

26.36X 

24.79* 

21.94* 

15.24* 

22.40* 

20.86* 

18.09* 

11.82* 

91.89* 

91.06* 

89.43* 

84.72X 

74.80* 

71.48* 

65.29X 

49.83* 

60.21* 

54.99* 

45 . 88* 

26.72* 

40.33* 

33.64* 

23.41* 

7.94* 

19.00* 

13.22* 

6.40* 

0.74* 

9.26* 

5.37* 

1.81* 

0.07X 

3.23* 

1.43* 

0.28* 

O.OOX 

1.94* 

0.75% 

0.11* 

0.00* 

91.45* 

90 . 47* 

88.53* 

83.01* 

73.02* 

69.18* 

62.10X 

45.00* 

57.38* 

51 .50* 

41.50* 

21.79* 

36.62* 

29.50* 

19.15X 

5.28* 

15.66* 

10.16* 

4.29* 

0.33* 

6.93* 

3.62* 

0.99* 

0.02* 

2.09* 

0.79* 

0.11* 

0.00* 

1.17* 

0.38* 

0.04* 

0.00* 

91.33* 

90.30* 

88.29* 

82.55* 

72.54* 

68.55* 

61.25* 

43.76* 

56.62* 

50.58* 

40.37* 

20.60* 

35.66* 

28.45* 

18.13* 

4.72X 

14.86* 

9.46* 

3.84* 

0.26* 

6.40* 

3.25* 

0.84* 

0.01* 

1.86* 

0.67* 

0.09* 

0.00* 

1.02* 

0.31* 

0.03X 

O.OOX 

91.27* 

90.23* 

88.18* 

82.33* 

72.31* 

68.27* 

60.85* 

43.19* 

56.27* 

50.15* 

39.85* 

20.07* 

35.22* 

27.98X 

17.66* 

4.48* 

14.49* 

9.14* 

3.65* 

0.23* 

6.16* 

3.09* 

0.78* 

0.01* 

1.76* 

0.62* 

0.08* 

0.00* 

0.95* 

0.28* 

0.03* 

0.00* 

86 
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TAB  LB  A.  2.b 


TRANSMITTANCE  IN  3.5  UM  BAND  -  MIDLATITUDE  SUMMER  PROFILE 

***  EMPIRICAL  *** 


Range  !  Visibility  in  Kiloneters 

in  Kb  !  40  30  20  15  10 


Wind=  0 
Wind=  0 
Wind=  0 
Wind=  0 
Wind=  0 
Wind=  0 
Mind=  0 
Wind=  0 

Mind=  5 
Wind=  5 
Wind=  5 
Wind=  5 
Wind=  5 
Wind=  5 
Wind=  5 
Wind=  5 

Wind=10 
lfind=  10 
Mind=10 
Wind=10 
Wind=10 
Wind=10 
Wind=10 
Nind=l0 

Wind=15 
Wind=15 
Wind=15 
Wind=15 
Wind=15 
Wind= 15 
Wind=15 
Wind=15 

Wind=20 

Wind=20 

Wind=20 

Wind=20 

Wind-20 

Wind=20 

Wind=20 

Wind=20 


0.2  : 
1.0  : 
2.0  : 

4.0  : 
8.0  : 
12.0  : 
i8.o  : 
21.0  : 

0.2  : 
l.o  : 
2.0  : 
4.0  : 
8.o  : 
12.0  : 
18.0  : 
21.0  : 

0.2  : 
l.o  : 
2.0  : 

4.0  : 

8.o  : 
12.0  : 
i8.o  : 
21.0  : 

0.2  : 
l.o  : 
2.0  : 

4.0  : 

8.o  : 
12.0  : 
18.0  : 
21.0  : 

0.2  : 
l.o  : 
2.0  : 

4.0  : 
8.0  : 
12.0  : 
18.0  : 
21.0  : 


86.48% 

82.52% 

77.82% 

69.21% 

54.74% 

43.30% 

30.46% 

25.55% 

85 . 70% 
78.85% 
71.06% 
57.70% 
38.05% 
25.10% 
13.44% 
9.84% 

85 . 49% 
77 . 90% 
69.36% 
54.97% 
34.54% 
21.70% 
10.81% 
7.63% 

85.44% 
77.66% 
68.92% 
54.29% 
33.68% 
20.90% 
10.21% 
7. 14% 

85.42% 

77.59% 

68.81% 

54.11% 

33.46% 

20.70% 

10.06% 

7.02% 


86.39% 

82.07% 

76.99% 

67.73% 

52.43% 

40.59% 

27.64% 

22.81% 

85.07% 
76.01% 
66.03% 
49.83% 
28.38% 
16. 16% 
6.95% 
4.55% 

84.79% 

74.78% 

63.90% 

46.67% 

24.90% 

13.28% 

5.17% 

3.23% 

84.72% 

74.46% 

63.36% 

45.89% 

24.06% 

12.62% 

4.79% 

2.95% 

84.70% 

74.38% 

63.23% 

45.68% 

23.85% 

12.45% 

4.70% 

2.89% 


86.35% 

81.88% 

76.62% 

67.09% 

51.45% 

39.45% 

26.49% 

21.70% 

84.11% 

71.83% 

58.97% 

39.74% 

18.05% 

8.20% 

2.51% 

1.39% 

83.73% 

70.20% 

56.32% 

36.25% 

15.02% 

6.22% 

1.66% 

0.86% 

83.63% 

69.78% 

55.65% 

35.40% 

14.32% 

5.79% 

1.49% 

0.76% 

83.60% 

69.68% 

55.48% 

35.18% 

14.15% 

5.69% 

1.45% 

0.73% 


86.26% 

81.47% 

75.86% 

65.77% 

49.44% 

37.16% 

24.22% 

19.55% 

83.02% 

67.29% 

51.74% 

30.60% 

10.70% 

3.74% 

0.77% 

0.35% 

82.52% 

65.26% 

48.67% 

27.08% 

8.38% 

2.59% 

0.45% 

0.19% 

82.39% 

64.75% 

47.91% 

26.23% 

7.86% 

2.36% 

0.39% 

0.16% 

82.35% 

64.61% 

47.71% 

26.02% 

7.74% 

2.30% 

0.37% 

0.15% 


86.09% 

80.65% 

74.34% 

63.16% 

45.60% 

32.92% 

20.19% 

15.81% 

80.69% 

58.34% 

38.89% 

17.29% 

3.42% 

0.67% 

0.06% 

0.02% 

79.93% 

55.65% 

35.39% 

14.32% 

2.34% 

0.38% 

0.03% 

0.01% 

79.73% 

54.97% 

34.54% 

13.63% 

2.12% 

0.33% 

0.02% 

0.01% 

79.68% 

54.80% 

34.32% 

13.46% 

2.07% 

0.32% 

0.02% 

0.00% 


5 


85.79% 

79.28% 

71.83% 

58.96% 

39.73% 

26.77% 

14.81% 

11.01% 

74.23% 

38.46% 

16.90% 

3.26% 

0.12% 

0.00% 

0.00% 

0.00% 

72.83% 

34.95% 

13.96% 

2.23% 

0.06% 

0.00% 

0.00% 

0.00% 

72.47% 

34.09% 

13.28% 

2.02% 

0.05% 

0.00% 

0.00% 

0.00% 

72.37% 

33.87% 

13.11% 

1.96% 

0.04% 

0.00% 

0.00% 

0.U0% 


n  *7 


-  TABLE  A. 3. a  - 

TRANSMITTANCE  IN  3.5  DM  BAND  -  MIDLATITDDB  WINTER  PROFILE 

***  LOWTRAN  DATA  *** 


Wind=5 


Wind=10 


Wind=15 


Wind=20 


Range  ! 


.2 


1 

2 
4 
8 

12.00 
18.00 
21.00 

0.2 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.2 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.2 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.2 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 


Visibility  in  Kiloaeters 


40 

30 

20 

15 

10 

5 

97.55* 

97.52* 

97.45* 

97.39* 

97.26* 

96.86* 

92 . 46* 

92.30* 

91.99* 

91.68* 

91.06* 

89.22* 

88.34* 

88.04* 

87.45* 

86.85X 

85.68* 

82.26* 

82.24* 

81.69* 

80.59* 

79.50* 

77.36* 

71.32* 

73.25* 

72.28* 

70.35* 

68.45* 

64.82* 

55.09* 

66.20* 

64 . 88* 

62.30* 

59.80* 

55.10* 

43.18* 

57.51* 

55.80* 

52.51* 

49.38* 

43.68* 

30.29X 

53.76* 

51.90* 

48.34* 

45.00* 

39.00* 

25.45* 

96.50* 

96.07* 

95.21* 

94.35* 

92.66* 

87.78* 

87.57* 

85 . 64* 

81.88* 

78.25* 

71.47* 

54.55* 

79.25* 

75.79* 

69.27* 

63.27* 

52.78* 

30.75* 

66. 19* 

60.54* 

50 . 58* 

42.18* 

29.36* 

9.96* 

47.45* 

39.69* 

27.70* 

19.27* 

9.34* 

1.08* 

34.51* 

26.40* 

15 . 40* 

8.93* 

3.01* 

0.12* 

21.65* 

14.49* 

6.45* 

2.85* 

0.56* 

0.00* 

17.19* 

10.76* 

4. 19* 

1.62* 

0.24* 

0.00* 

96.28* 

95.77* 

94.75* 

93.73* 

91 . 72* 

85.99X 

86.58* 

84.31* 

79.90* 

75 . 69* 

67 . 94* 

49.21* 

77.47* 

73.45* 

65 . 98* 

59.21* 

47.69* 

25.02* 

63.24* 

56.86* 

45 . 87* 

36.94* 

23.97* 

6.60* 

43.32* 

35.01* 

22.79* 

14.78* 

6.22X 

0.47* 

30.11* 

21.87* 

11.49* 

6.00* 

1.64X 

0.03* 

17.64* 

10.92* 

4.16X 

1.57* 

0.22* 

0.00* 

13.54* 

7.74* 

2.51* 

0.81* 

0.08* 

0.00* 

96.22* 

95.69* 

94.62* 

93.56* 

91.47* 

85.51* 

86.31* 

83.95* 

79.37* 

75.01* 

67.00* 

47.83* 

76.98* 

72.82* 

65.10* 

58.13* 

46.38* 

23.64* 

62.46* 

55.88* 

44.66* 

35.62* 

22.67* 

5.89* 

42.25* 

33.82* 

21.60* 

13.74* 

5.57* 

0.38* 

29.00* 

20.77* 

10.60* 

5.38* 

1.39* 

0.02* 

16.68* 

10.11* 

3.69* 

1.33* 

0.17* 

0.00* 

12.68* 

7.07* 

2.18* 

0.66* 

0.06* 

O.OOX 

96.19* 

95.65* 

94.56* 

93.48* 

91.35* 

85 . 28* 

86. 19* 

83.78* 

79. 12* 

74.69* 

66.56* 

47.20* 

76.76* 

72.52* 

64.69* 

57.64* 

45.78* 

23.02* 

62.09* 

55.43* 

44.10* 

35.01* 

22.09* 

5.59* 

41.76* 

33.28* 

21.07* 

13.28X 

5.28* 

0.34* 

28.49* 

20.27* 

10.21* 

5.11* 

1.28* 

0.02* 

16.24* 

9.75* 

3.48* 

1.23* 

0.16* 

0.00* 

12.30* 

6.78* 

2.04* 

0.61* 

0.05* 

0.00* 

TAB LB  4. 3.b 


TRANSMITTANCE  IN  3.5  (JM  BAND  -  MIDLATITUDE  WINTER  PROFILE 

***  EMPIRICAL  *** 


V 

Range  ! 

Visibility  in 

Kiloaeters 

{ 

in  Ka  ! 

40 

30 

20 

15 

10 

5 

\ 

Wind=  0 

0.2  : 

93.30* 

93.24* 

93.23* 

93. 15* 

92.97* 

92.65X 

-! 

Wind=  0 

1.0  : 

91.21* 

90.93* 

90.87* 

90.47* 

89.59* 

88.08* 

Wind=  0 

2.0  : 

88.66* 

88. 12* 

87 . 99* 

87.22* 

85.54* 

82.68* 

Wind=  0 

4.0  : 

83.77* 

82.75* 

82.52* 

81.08* 

77 . 99* 

72.85* 

Wind=  0 

8.o  : 

74.78* 

72.98* 

72.57* 

70.06* 

64.82* 

56.56* 

I 

Wind=  0 

12.0  : 

66.76* 

64.37* 

63.83* 

60 . 54* 

53.87* 

43.92* 

1 

Wind=  0 

18.0  : 

56.32* 

53.31* 

52.64* 

48.63* 

40.82* 

30.05* 

•; 

Wind=  0 

21.0  : 

51.72* 

48.52* 

47.81* 

43.59* 

35.53* 

24.85* 

Wind=  5 

0.2  : 

92.32* 

91.65* 

90.61* 

89.09* 

85.29* 

73.25* 

>.* 

Wind=  5 

l.o  : 

86.51* 

83 . 44* 

78.80* 

72.43* 

58.24* 

27.21* 

h 

Wind=  5 

2.0  : 

79.75* 

74.19* 

66. 18* 

55.90* 

36.15* 

7.89* 

$ 

Wind=  5 

4.0  : 

67 . 79* 

58.67* 

46.67* 

33.31* 

13.93* 

0.66* 

Wind=  5 

8.o  : 

48.97* 

36.68* 

23.22X 

11.82* 

2.07* 

O.OOX 

Wind=  5 

12.0  ; 

35.38* 

22.94* 

11.55* 

4.20X 

0.31* 

o.oox 

Wind=  5 

i8.o  : 

21.73* 

11.34* 

4.05* 

0.89* 

0.02* 

0.00* 

> 

Wind=  5 

21.0  : 

17.03* 

7.97* 

2.40* 

0.41* 

0.00* 

O.OOX 

i 

Wind=10 

0.2  ; 

92.09* 

91.33* 

90.14* 

88.41* 

84.08* 

70.58X 

Wind=10 

l.o  : 

85 . 46* 

81.98* 

76.77* 

69.67* 

54.22* 

22.60X 

„  • 

Wind=10 

2.0  i 

77.83* 

71.63* 

62.81* 

51.73* 

31.34* 

5.44* 

V’ 

Wind=10 

4.0  : 

64.57* 

54.68* 

42.04* 

28.52* 

10.46* 

0.32* 

Wind=10 

8.o  ; 

44.43* 

31.87* 

18.83* 

8.67* 

1.17* 

0.00* 

■ 

Wind=10 

12.0  : 

30.57X 

18.57* 

8.44* 

2.63X 

0.13* 

0.00* 

1 

Wind=10 

i8.o  : 

17.45* 

8.26X 

2.53* 

0.44* 

0.00* 

0.00* 

w 

Wind=10 

21.0  : 

13. 18* 

5.51* 

1.39* 

0.  18* 

0.00* 

0.00* 

Wind=15 

0.2  : 

92.04* 

91.26* 

90.03* 

88.25* 

83.81* 

69.99* 

Wind= 15 

l.o  : 

85.22* 

81.65* 

76.31* 

69.05* 

53.34X 

21.66* 

Wind=15 

2.0  : 

77.40* 

71.06* 

62.06* 

50.82* 

30.33* 

5.00* 

I 

Wind= 15 

4.0  : 

63.85* 

53.81* 

41.04* 

27.53* 

9.80* 

0.27* 

», 

Wind=15 

8.0  : 

43.45* 

30.86* 

17.95* 

8.07* 

1.02* 

0.00* 

k  • 

Wind=15 

12.0  : 

29.57* 

17.70* 

7.85* 

2.37* 

0.11* 

0.00* 

Wind=15 

18.0  : 

16.60* 

7.69* 

2.27X 

0.38* 

0.00* 

O.OOX 

Wind=15 

21.0  : 

12.43* 

5.06X 

1.22* 

0.15* 

0.00* 

O.OOX 

«-  - 

B 

Wind=20 

0.2  : 

92.03* 

91.24* 

90.01* 

88.21* 

83.75* 

69.85* 

Wind=20 

l.o  : 

85.17* 

81.58X 

76.20X 

68.91* 

53. 15* 

21.45* 

•* 

Wind=20 

2.0  ! 

77.30* 

70.92* 

61.89* 

50.61* 

30. 10X 

4.91* 

»  * 

Wind=20 

4.0  : 

63.69X 

53.61* 

40.82* 

27.30* 

9.66* 

0.26* 

►  * 

Wind=20 

8.o  : 

43.23* 

30.63* 

17.76* 

7.94* 

0.99* 

0.00* 

A 

Wind=20 

12.0  : 

29.34* 

17.50* 

7.73* 

2.31* 

0. 10* 

0.00* 

K 

Wind=20 

i8.o  : 

16.41* 

7.56* 

2.22* 

0.36* 

0.00* 

0.00* 

•J. 

Wind=20 

21.0  : 

12.27* 

4.97* 

1.19* 

0.14* 

0.00* 

0.00* 

» . 

»* 

V 

\ 

89 
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-  TABLE  A. 4. a  - 

TRANSMITTANCE  IN  3.5  UN  BAND  -  SUB-ARCTIC  SUMMER  PROFILE 

***  LOWTRAN  DATA  tt* 


Range  ! 
in  Ka  ! 

40  Ka 

Visibility  in  Kiloaeters 

30  Kb  20  Kb  15  Ka  10  Kb 

5  Kb 

Wind=0 

0.20 

95.52X 

95.49k 

95.43k 

95.36k 

95.23k 

94 . 85k 

1. 00 

87 . 64k 

87.49k 

87.20k 

86.91k 

86.32k 

84.60k 

2.00 

81.41k 

81.14k 

80.60k 

80.05k 

78.98k 

75 . 86k 

4.00 

72.23k 

71.75k 

70.79k 

69.84k 

67 . 98k 

62.72k 

8.00 

59.00k 

58.22k 

56.68k 

55.17k 

52.26k 

44.48k 

12.00 

49.49k 

48.51k 

46.60k 

44.74k 

41.26k 

32.40k 

18.00 

38.90k 

37.75k 

35.54k 

33.44k 

29.61k 

20 . 60k 

21.00 

34.70k 

33.50k 

31.23k 

29.09k 

25.24k 

16.53k 

Wind=5 

0.20 

94.49k 

94.07k 

93.23k 

92.38k 

90.73k 

85.95k 

1.00 

83.00k 

81.17k 

77.60k 

74.16k 

67.73k 

51.69k 

2.00 

73.02k 

69.84k 

63.83k 

58.29k 

48.63k 

28.32k 

4.00 

58.11k 

53.15k 

44.40k 

37.03k 

25.77k 

8.74k 

8.00 

38.19k 

31.95k 

22.30k 

15.51k 

7.51k 

0.86k 

12.00 

25.78k 

19.72k 

11.50k 

6.67k 

2.25k 

0.09k 

18.00 

14.62k 

9.79k 

4.36k 

1.92k 

0.38k 

0.00k 

21.00 

11.08k 

6.93k 

2.70k 

1.04k 

0. 16k 

0.00k 

Wind=10  0.20 

94.28k 

93.78k 

92.78k 

91.79k 

89.83k 

84.23k 

1.00 

82 . 08k 

79.93k 

75.77k 

71.79k 

64.45k 

46.73k 

2.00 

71.42k 

67.72k 

60.85k 

54.63k 

44.03k 

23.15k 

4.00 

55.58k 

49.98k 

40.36k 

32.52k 

21.13k 

5.84k 

8.00 

34.94k 

28.26k 

18.42k 

11.96k 

5.05k 

0.39k 

12.00 

22.55k 

16.40k 

8.63k 

4.52k 

1.24k 

0.03k 

18.00 

11.97k 

7.42k 

2.83k 

1.07k 

0.15k 

0.00k 

21.00 

8.77k 

5.02k 

1.63k 

0.53k 

0.05k 

o.ook 

Wind=15  0.20 

94.22k 

93.70k 

92.66k 

91.62k 

89.58k 

83.77k 

1.00 

81.83k 

79.60k 

75.28k 

71.15k 

63.58k 

45.46k 

2.00 

70.98k 

67.16k 

60.06k 

53.66k 

42 . 85k 

21.90k 

4.00 

54.91k 

49.15k 

39.31k 

31.38k 

20.01k 

5.23k 

8.00 

34.10k 

27.32k 

17.48k 

11.14k 

4.53k 

0.31k 

12.00 

21.74k 

15.59k 

7.98k 

4.06k 

1.05k 

0.02k 

18.00 

11.33k 

6.88k 

2.52k 

0.91k 

0.12k 

o.ook 

21.00 

8.23k 

4.60k 

1.42k 

0.44k 

0.04k 

o.ook 

Wind=20  0.20 

94.20k 

93.67k 

92.60k 

91.55k 

89.47k 

83.55k 

1.00 

81.72k 

79.44k 

75.05k 

70.86k 

63.18k 

44.87k 

2.00 

70.78k 

66.89k 

59.70k 

53.22k 

42.31k 

21.34k 

4.00 

54.60k 

48.76k 

38.84k 

30.86k 

19.51k 

4.96k 

8.00 

33.71k 

26.89k 

17.06k 

10.77k 

4.30k 

0.28k 

12.00 

21.37k 

15.23k 

7.69k 

3.86k 

0.97k 

0.02k 

18.00 

11.04k 

6.64k 

2.38k 

0.85k 

0.11k 

o.ook 

21.00 

7.98k 

4.41k 

1.34k 

0.40k 

0.04k 

o.ook 

90 


-  TABLE  A.4.b  - 

TRANSMITTANCE  IN  3.5  DM  BAND  -  SUB-ARCTIC  SUMMER  PROFILE 

***  EMPIRICAL  *** 


Range  ! 
in  Ka  !  40 


Visibility  in  Kiloaeters 
30  20  15  10 


92 . 92% 

92.87% 

92.85% 

92.77% 

92.59% 

92.28% 

89.38% 

89.11% 

89.05% 

88.66% 

87.81% 

86.34% 

85.13% 

84.63% 

84.51% 

83.77% 

82.17% 

79.44% 

77.24% 

76.33% 

76.11% 

74.79% 

71.96% 

67 . 26% 

63.59% 

62.09% 

61.74% 

59.62% 

55.19% 

48.22% 

52.35% 

50.51% 

50.08% 

47.52% 

42.32% 

34.57% 

39.10% 

37.06% 

36.59% 

33.82% 

28.43% 

20.98% 

33.79% 

31.75% 

31.28% 

28.53% 

23.30% 

16.35% 

91.95% 

91.33% 

90.33% 

89.08% 

86.25% 

78.02% 

84.81% 

81.96% 

77.59% 

72.38% 

61.59% 

37  30% 

76.65% 

71.59% 

64.16% 

55.83% 

40.43% 

14.83% 

62.62% 

54.63% 

43.88% 

33.22% 

17.42% 

2.34% 

41.79% 

31.80% 

20.52% 

11.76% 

3.23% 

0.06% 

27.89% 

18.52% 

9.60% 

4. 17% 

0.60% 

0.00% 

15.20% 

8.23% 

3.07% 

0.88% 

0.05% 

0.00% 

11.23% 

5.48% 

1.74% 

0.40% 

0.01% 

0.00% 

91.73% 

91.04% 

89.93% 

88.54% 

85.40% 

76.34% 

83 . 80% 

80 . 67% 

75.87% 

70.20% 

58.61% 

33.46% 

74.85% 

69.35% 

61.36% 

52.52% 

36.61% 

11.93% 

59.70% 

51.25% 

40.12% 

29.40% 

14.29% 

1.52% 

37.99% 

28.00% 

17.15% 

9.21% 

2.17% 

0.02% 

24. 17% 

15.29% 

7.34% 

2.89% 

0.33% 

0.00% 

12.27% 

6.17% 

2.05% 

0.51% 

0.02% 

0.00% 

8.74% 

3.92% 

1.08% 

0.21% 

0.00% 

0.00% 

91.68% 

90.97% 

89.84% 

88.42% 

85.21% 

75.97% 

83.58% 

80.38% 

75.50% 

69.72% 

57.97% 

32 . 66% 

74.45% 

68.86% 

60.75% 

51.81% 

35.81% 

11.37% 

59.07% 

50.53% 

39.33% 

28.61% 

13.67% 

1.38% 

37.19% 

27.21% 

16.48% 

8.72% 

1.99% 

0.02% 

23.41% 

14.66% 

6.91% 

2.66% 

0.29% 

0.00% 

11.70% 

5.79% 

1.87% 

0.45% 

0.02% 

0.00% 

8.27% 

3.64% 

0.98% 

0. 18% 

0.00% 

0.00% 

91.67% 

90 . 96% 

89.82% 

88.39% 

85.17% 

75.89% 

83.53% 

80.32% 

75.41% 

69.62% 

57.82% 

32.48% 

74.36% 

68.75% 

60.61% 

51.65% 

35.63% 

11.24% 

58.93% 

50.37% 

39. 15% 

28.43% 

13.53% 

1.35% 

37.01% 

27.04% 

16.34% 

8.62% 

1.95% 

0.02% 

23.25% 

14.52% 

6.82% 

2.61% 

0.28% 

0.00% 

11.57% 

5.71% 

1.84% 

0.44% 

0.02% 

0.00% 

8. 16% 

3.58% 

0.95% 

0.18% 

0.00% 

0.00% 

TABLE  A. 5. a 


TRANSMITTANCE  IN  3.5  UM  BAND  -  SUB-ARCTIC  WINTER 

***  LOWTRAN  DATA  *** 


Wind=0 


Wind=5 


Wind=10 


Wind=15 


Wind=20 


Visibility  in  Eiloaeters 
40  30  20  15  10 


Range  ! 
in  la  ! 


0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 


98.73* 

95.78* 

93.00* 

88.76* 

82.56* 

77.53* 

71.13* 

68.30* 

97.67* 

90.73* 

83.45* 

71.47* 

53.53* 

40.47* 

26.83* 

21.90* 

97.44* 
89.67* 
81.50* 
68. 18X 
48.70* 
35.13* 
21.70* 
17.09* 

97 . 38* 
89.37* 
80.97* 
67.29* 
47.44* 
33.78* 
20.45* 
15.96* 

97.35* 

89.24* 

80.72* 

66.88* 

46.86* 

33.16* 

19.90* 

15.45* 


98.70* 
95.62* 
92.68* 
88. 16* 
81.44* 
75.96* 
68.98* 
65.90* 

97.24* 
88.73* 
79.81* 
65.38* 
44.78* 
30 . 98* 
17.96* 
13.71* 

96.92* 

87.30* 

77.26* 

61.25* 

39.32* 

25.48* 

13.40* 

9.74* 

96.83* 
86.90* 
76.56* 
60.16* 
37.92* 
24. 13* 
12.35* 
8.86* 

96.79* 

86.72* 

76.24* 

59.65* 

37.28* 

23.53* 

11.89* 

8.47* 


98.63* 
95.29* 
92.05* 
86.96* 
79.23* 
72.89* 
64 . 84* 
61.31* 

96.37* 
84.83* 
72 . 96* 
54.63* 
31.27* 
18.07* 
8.00X 
5.34* 

95.88* 
82.71* 
69.36* 
49.37* 
25.54* 
13.34* 
5.08* 
3. 14X 

95.74X 

82.13* 

68.39* 

48.00* 

24.14* 

12.26* 

4.47* 

2.71* 

95 . 68* 
81.86* 
67.94* 
47.37* 
23.51* 
11.78* 
4.21* 
2.53* 


98.56* 

94.96* 

91.41* 

85.76* 

77.06* 

69.92* 

60.92* 

57.01* 

95.50* 

81.08* 

66.64* 

45.57* 

21.76* 

10.49* 

3.54* 

2.06* 

94 . 84* 
78.33* 
62.20* 
39.71* 
16.52* 
6.94* 
1.91* 
1.00* 

94 . 66* 
77.59* 
61.02* 
38.22* 
15.30* 
6.  19* 
1.60* 
0.82* 

94.58* 

77.24* 

60.48* 

37.54* 

14.76* 

5.86* 

1.48X 

0.74* 


98.43* 
94.30* 
90. 15* 
83.42* 
72.91* 
64.34* 
53.78* 
49.29* 

93.78* 
74.06* 
55 . 60* 
31.73* 
10.55* 
3.54* 
0.69X 
0.31* 

92.80* 
70 . 26* 
50.05X 
25.71* 
6.92* 
1.88* 
0.27* 
0.10* 

92.53* 
69.24* 
48.61* 
24.25* 
6. 16* 
1.58* 
0.21* 
0.08* 

92.40* 

68.77* 

47.94* 

23.59* 

5.83* 

1.46* 

0.18* 

0.07* 


PROFILE 


5 


98.02* 

92.38* 

86.51* 

76.80* 

61.81* 

50.22* 

37.08* 

31.94* 

88.85* 

56.54* 

32.41* 

10.78* 

1.22* 

0.14* 

0.01* 

0.00* 

86.97* 

50.81* 

26.17* 

7.03* 

0.52* 

0.04* 

0.00* 

0.00* 

86.46* 

49.32* 

24.66* 

6.24* 

0.41* 

0.03* 

0.00* 

0.00* 

86.22* 

48.64* 

23.98* 

5.90* 

0.36* 

0.02* 

0.00* 

0.00* 
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TABLE  A . 5 . b 


TRANSMITTANCE  IN  3.S  UM  BAND  -  SUBARCTIC  WINTER  PROFILE 

***  EMPIRICAL  «** 


Range  !  Visibility  in  Kiloneters 


in  Ka  ! 

40 

30 

20 

15 

10 

5 

Wind= 

0 

0.2 

96.02X 

95 . 96% 

95 . 95% 

95.86% 

95 . 68% 

95.34% 

Wind= 

0 

1.0 

94.69% 

94.40% 

94.33% 

93.91% 

92 . 98% 

91.39% 

Wind= 

0 

2.0 

93.05% 

92.47% 

92.34% 

91.52% 

89.73% 

86.67% 

Wind= 

0 

4.0 

89.85% 

88.74% 

88.49% 

86.92% 

83.55% 

77.96% 

Wind= 

0 

8.0 

83.78% 

81.73% 

81.26% 

78.41% 

72.44% 

63.08% 

Wind= 

0 

12.0 

78.12% 

75.27% 

74.63% 

70.73% 

62.81% 

51.03% 

Wind= 

0 

18.0 

70.34% 

66.52% 

65.67% 

60 . 59% 

50.71% 

37.14% 

Wind= 

0 

21.0 

66.74% 

62.54% 

61.61% 

56.09% 

45.57% 

31.68% 

Wind= 

5 

0.2 

95 . 02% 

94.28% 

93.17% 

91.28% 

86.18% 

69.38% 

Wind= 

5 

1.0 

89.86% 

86.39% 

81.42% 

73.51% 

55. 13% 

18.65% 

Wind= 

5 

2.0 

83.80% 

77.46% 

68.80% 

56.08% 

31.54% 

3.61% 

Wind= 

5 

4.0 

72.87% 

62.26% 

49. 12% 

32.64% 

10.32% 

0.14% 

Wind= 

5 

8.0 

55.11% 

40.23% 

25.04% 

11.06% 

1.11% 

0.00% 

Wind= 

5 

12.0 

41.68% 

26.00% 

12.76% 

3.75% 

0.12% 

0.00% 

Wind= 

5 

18.0 

27.41% 

13.50% 

4.65% 

0.74% 

0.00% 

0.00% 

Wind= 

5 

21.0 

22.23% 

9.73% 

2.80% 

0.33% 

0.00% 

0.00% 

Wind= 

10 

0.2 

94.78% 

93.90% 

92.60% 

90 . 38% 

84.44% 

65.34% 

Wind  = 

10 

1.0 

88.72% 

84.69% 

78.95% 

69 . 96% 

49.78% 

13.81% 

Wind= 

10 

2.0 

81.68% 

74.43% 

64.69% 

50.80% 

25.71% 

1.98% 

Wind  = 

10 

4.0 

69.24% 

57.49% 

43.43% 

26.78% 

6.86% 

0.04% 

Wind= 

10 

8.0 

49.76% 

34.30% 

19.57% 

7.44% 

0.49% 

0.00% 

Wind= 

10 

12.0 

35.76% 

20 . 46% 

8.82% 

2.07% 

0.03% 

0.00% 

Wind= 

10 

18.0 

21.78% 

9.43% 

2.67% 

0.30% 

0.00% 

0.00% 

Wind= 

10 

21.0 

17.00% 

6.40% 

1.47% 

0.12% 

0.00% 

0.00% 

Wind= 

15 

0.2 

94.72% 

93.81% 

92.46% 

90.17% 

84.02% 

64.39% 

Wind= 

15 

1.0 

88.44% 

84.28% 

78.36% 

69. 13% 

48.56% 

12.84% 

Wind= 

15 

2.0 

81.18% 

73.71% 

63.73% 

49.59% 

24.47% 

1.71% 

Wind= 

15 

4.0 

68.39% 

56.38% 

42.14% 

25.52% 

6.21% 

0.03% 

Wind= 

15 

8.0 

48.54% 

32.99% 

18.43% 

6.76% 

0.40% 

0.00% 

Wind= 

15 

12.0 

34.45% 

19.30% 

8.06% 

1.79% 

0.03% 

0.00% 

Wind= 

15 

18.0 

20.60% 

8.64% 

2.33% 

0.24% 

0.00% 

0.00% 

Wind= 

15 

21.0 

15.93% 

5.78% 

1.25% 

0.09% 

0.00% 

0.00% 

Wind= 

20 

0.2 

94.71% 

93.79% 

92.42% 

90.11% 

83.91% 

64.16% 

Wind= 

20 

1.0 

88.38% 

84.18% 

78.22% 

68.92% 

48.26% 

12.61% 

Wind= 

20 

2.0 

81.06% 

73.53% 

63.49% 

49.30% 

24. 17% 

1.65% 

Wind= 

20 

4.0 

68.18% 

56.11% 

41.84% 

25.22% 

6.06% 

0.03% 

Wind= 

20 

8.0 

48.24% 

32.68% 

18.17% 

6.60% 

0.38% 

0.00% 

Wind= 

20 

12.0 

34.14% 

19.03% 

7.89% 

1.73% 

0.02% 

0.00% 

Wind= 

20 

18.0 

20.32% 

8.46% 

2.26% 

0.23% 

0.00% 

0.00% 

Wind= 

20 

21.0 

15.67% 

5.64% 

1.21% 

0.08% 

0.00% 

0.00% 
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TABLE  A.  6.  a 


TRANSMITTANCE  IN  3.5  UM  BAND  -  STANDARD  ATMOSPHBRE  PROFILE 

tt*  LOWTRAN  DATA  *#* 

Range  !  Visibility  in  Kiloneters 

in  Kn  !  40  30  20  15  10  5 


Wind=0  0.20  : 

1.00  : 
2.00  : 
4.00  : 
8.oo  : 
12.00  : 
18.00  : 
21.00  : 

Wind=5  0.20  ! 

1.00  : 
2.00  : 
4.00  : 
8.oo  : 
12.00  : 
i8.oo  : 
21.00  : 

Wind= 10  0.20  ! 

1.00  ! 
2.oo  : 
4.00  : 
8.oo  : 
12.00  : 
I8.00  : 
21.00  : 

Wind=15  0.20  ! 

1.00  : 

2.00  : 

4.oo  : 

8.oo  : 
12.00  : 
18.00  : 
21.00  : 

Wind=20  0.20  : 

1.00  ! 
2.00  : 
4.00  : 
8.oo  : 
12.00  : 
i8.oo  : 
21.00  : 


96.67%  96.64% 
90.27%  90.13% 
85.25%  84.99% 
77.78%  77.32% 
67.02%  66.22% 
58.64%  57.60% 
49.00%  47.70% 
45.01%  43.62% 

95.61%  95.18% 
85.43%  83.53% 
76.35%  73.00% 
62.39%  57.04% 
43.12%  36.03% 
30.26%  23.12% 
18.17%  12.13% 
14.14%  8.83% 

95 . 45%  94 . 96% 
84.72%  82.57% 
75.08%  71.33% 
60.34%  54.46% 
40.33%  32.86% 
27.37%  20.13% 
15.63%  9.86% 
11.87%  6.93% 

95.40%  94.90% 
84.53%  82.33% 
74.76%  70.91% 
59.82%  53.82% 
39.64%  32.08% 
26.67%  19.42% 
15.03%  9.34% 
11.34%  6.51% 

95.39%  94.88% 
84.45%  82.21% 
74.61%  70.71% 
59.58%  53.52% 
39.32%  31.73% 
26.36%  19.10% 
14.77%  9.11% 
11.10%  6.32% 


96.58% 

96.52% 

89.86% 

89.59% 

84.48% 

83.97% 

76.40% 

75.47% 

64.65% 

63.09% 

55.55% 

53.57% 

45.19% 

42.78% 

40.95% 

38.42% 

94.32% 

93.46% 

79.84% 

76.28% 

66.68% 

60.87% 

47.60% 

39.66% 

25.10% 

17.42% 

13.44% 

7.77% 

5.37% 

2.36% 

3.42% 

1.31% 

93.98% 

93.01% 

78.42% 

74.44% 

64 . 34% 

57.97% 

44.30% 

35.97% 

21.74% 

14.33% 

10.84% 

5.80% 

3.89% 

1.52% 

2.34% 

0.79% 

93.90% 

92 . 89% 

78.06% 

73.97% 

63.74% 

57.24% 

43.49% 

35.07% 

20.95% 

13.63% 

10.25% 

5.38% 

3.58% 

1.36% 

2.13% 

0.69% 

93.86% 

92.84% 

77.89% 

73.76% 

63.47% 

56.91% 

43.12% 

34.67% 

20.59% 

13.31% 

9.99% 

5.19% 

3.45% 

1.29% 

2.03% 

0.65% 

96.40% 

96.05% 

89.05% 

87.45% 

82.96% 

80.01% 

73.67% 

68.52% 

60.11% 

52.01% 

49.81% 

40.09% 

38.36% 

27.70% 

33.82% 

23.13% 

91.78% 

86.92% 

69.63% 

53.06% 

50.73% 

24.95% 

27.54% 

9.29% 

8.40% 

0.96% 

2.60% 

0.10% 

0.46% 

0.00% 

0.19% 

0.00% 

91.09% 

85.61% 

67 . 09% 

49.19% 

47.08% 

25.31% 

23.73% 

6.86% 

6.24% 

0.52% 

1.66% 

0.04% 

0.23% 

0.00% 

0.09% 

0.00% 

90.92% 

85.28% 

66.44% 

48.23% 

46. 18% 

24.34% 

22.83% 

6.34% 

5.77% 

0.45% 

1.48% 

0.03% 

0.20% 

0.00% 

0.07% 

0.00% 

90 . 84% 

85.12% 

66. 15% 

47 . 80% 

45.77% 

23.90% 

22.43% 

6. 12% 

5.57% 

0.41% 

1.41% 

0.03% 

0.18% 

0.00% 

0.07% 

0.00% 
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TABLE  A. 6.b 


TRANSMITTANCE  IN  3.5  DM  BAND  -  STANDARD  ATMOSPHERE 

***  EMPIRICAL  tt* 

Range  !  Visibility  in  Eiloaeters 


in  !■  !  40  30 


*find= 

0 

0.2 

91.41* 

91.36* 

Wind= 

0 

1.0 

88.78* 

88.53* 

Wind= 

0 

2.0 

85 . 59* 

85.13* 

Wind= 

0 

4.0 

79.56* 

78.70* 

Wind= 

0 

8.0 

68.73* 

67.26* 

Wind= 

0 

12.0 

59.38* 

57.48* 

Wind= 

0 

18.0 

47.69* 

45.42* 

Wind= 

0 

21.0 

42.74* 

40.37* 

Wind= 

5 

0.2 

90.45* 

89.84* 

Wind= 

5 

1.0 

84.21* 

81.41* 

Wind= 

5 

2.0 

77.01* 

71.98* 

Hind= 

5 

4.0 

64.41* 

56.26* 

Wind= 

5 

8.0 

45.06* 

34.38* 

Wind= 

5 

12.0 

31.52* 

21.01* 

Wind= 

5 

18.0 

18.44* 

10.03* 

Wind= 

5 

21.0 

14.11* 

6.93* 

Wind= 

10 

0.2 

90.27* 

89.60* 

Wind  = 

10 

1.0 

83.39* 

80.33* 

Wind= 

10 

2.0 

75.53* 

70.07* 

Wind= 

10 

4.0 

61.95* 

53.32* 

Wind= 

10 

8.0 

41.68* 

30.88* 

Wind= 

10 

12.0 

28.04* 

17.88* 

Wind= 

10 

18.0 

15.48* 

7.88* 

Wind= 

10 

21.0 

11.50* 

5.23* 

Wind= 

15 

0.2 

90.24X 

89.55X 

Wind  = 

15 

1.0 

83.25* 

80.13* 

Wind= 

15 

2.0 

75.26* 

69.73* 

Wind= 

15 

4.0 

61.51* 

52.80* 

Wind= 

15 

8.0 

41.09* 

30.28* 

Wind= 

15 

12.0 

27.45* 

17.36X 

Wind= 

15 

18.0 

14.99* 

7.54X 

Wind= 

15 

21.0 

11.08* 

4.97* 

Wind= 

20 

0.2 

90.24* 

89.55* 

Wind= 

20 

1.0 

83.22* 

80.09* 

Wind= 

20 

2.0 

75.21* 

69 . 67* 

Wind= 

20 

4.0 

61.43* 

52.71* 

Nind= 

20 

8.0 

40.99* 

30. 17* 

Wind  = 

20 

12.0 

27.34* 

17.27* 

Wind= 

20 

18.0 

14.90* 

7.48* 

Wind= 

20 

21.0 

11.00* 

4.92* 

20 

15 

10 

91.35* 

91.28* 

91.12* 

88.48* 

88. 13* 

87.37* 

85.02* 

84.36* 

82.91X 

78.50* 

77.28* 

74.65* 

66.92* 

64 . 86* 

60.52* 

57.06* 

54.44* 

49.06* 

44.91* 

41.86* 

35.81* 

39.85* 

36.70* 

30.60* 

88.86* 

87.58* 

84.55* 

77.08* 

71.67* 

60.11* 

64.53* 

55.79* 

39.23* 

45.22* 

33.80* 

16.72* 

22.20X 

12.40* 

3.04* 

10.90* 

4.55* 

0.55* 

3.75* 

1.01* 

0.04* 

2.20* 

0.48* 

0.01* 

88.52* 

87. 10X 

83.77X 

75.60* 

69.74* 

57.37* 

62.07* 

52.82* 

35.75* 

41.84* 

30.30* 

13.88X 

19.01* 

9.97X 

2.09* 

8.64* 

3.28* 

0.32* 

2.65X 

0.62* 

0.02* 

1.46* 

0.27* 

0.00* 

88.46* 

87.01* 

83.63* 

75.33X 

69.39* 

56.89* 

61.63* 

52.29* 

35.15* 

41.25* 

29.70X 

13.42* 

18.48X 

9.58* 

1.95* 

8.28* 

3.09* 

0.28* 

2.48* 

0.57* 

0.025* 

1.36* 

0.24* 

0.00* 

88.45* 

87.00* 

83.60* 

75.28X 

69.33* 

56.80* 

61.55* 

52.20* 

35.01* 

41.15* 

29.59* 

13.33* 

18.39* 

9.51* 

1.93* 

8.22* 

3.06* 

0.28* 

2.45* 

0.56* 

0.02* 

1.34* 

0.24* 

0.00* 

PROFILE 
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90 . 84* 
86.06* 
80.43* 
70.26* 
53.61* 
40.90* 
27.26* 
22.26* 

75.35* 

33.78* 

12.39* 

1.67* 

0.03* 

0.00* 

0.00* 

0.00* 

73.72* 
30.28* 
9. 96* 
1.08* 
0.01* 
0.00* 
0.00* 
0.00* 

73.42* 

29.68* 

9.57* 

0.99* 

0.01* 

0.00* 

0.00* 

0.00* 

73.37* 

29.57* 

9.50* 

0.98* 

0.01* 

0.00* 

0.00* 

0.00* 
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TABLE  A. 7 


TABLE  OF  A— COEFFICIENTS  FOR  SIX  LOWTRAN  PROFILES  AT  3.5  urn. 


Tropical  Profile 


Visibi 1 ity 

40 

30 

20 

15 

lO 

Mind  =  O 

0.856781 

0.856471 

0.856390 

0.856255 

0.856115 

Wind  -  5 

0.856363 

0.856405 

0.856220 

0.857041 

0.859499 

Wind  =  10 

0.855773 

0. 856587 

0.856102 

0.855903 

0.851515 

Wind  =  15 

0.856526 

0.856341 

0.855761 

0.856530 

0.862983 

Wind  =  20 

0.856982 

0.857127 

0.856171 

0.857050 

0.851565 

Midlatitude  Summer  Profile 

Vi sibi 1 ity 

40 

30 

20 

15 

10 

Wind  =  O 

0.875788 

0. 875808 

0.875754 

0.B75834 

0.875812 

Wind  -  5 

0.875747 

0.875824 

0.875699 

0.875477 

0.877011 

Wind  ■  lO 

0.875820 

0.875716 

0.875650 

0.874656 

0.874052 

Wind  -  15 

O. 875753 

0.875685 

0.875366 

0.876011 

0.872420 

Wind  =  20 

0.875777 

0.875885 

0.875886 

0.B77431 

0.859462 

Midlatitude  Winter  Profile 

Visibi 1 ity 

40 

30 

20 

15 

10 

0.938293 

0.938300 

0.938272 

0.938302 

0.938302 

0.938338 

0.938338 

0.938284 

0.937979 

0.938643 

0.938264 

0.938385 

0.938206 

0.937628 

0.942839 

0.938298 

0.938291 

0.938232 

0.939094 

O. 942642 

0.938236 

0.938222 

0.938328 

0.937932 

0.945244 
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TABLE  OF  A— COEFFICIENTS  FOR  SIX  LOWTRAN  PROFILES  AT  3.5  um. 


Subarctic  Summer  Profile 


Visibi 1 ity 

40 

30 

20 

15 

io 

Wind  =  0 

0. 900942 

0.900967 

0.900928 

0.900916 

0.900932 

Wind  =  5 

0.900908 

0.900991 

0.900808 

0.901004 

0.897534 

Wind  =  10 

0.900920 

0.900967 

0.901247 

0.900275 

0.910749 

Wind  =  is 

0.900901 

0.900910 

0.901037 

0.900193 

0.903370 

Wind  =  20 

0.901048 

0.900953 

0.900717 

0.900709 

0.889788 

Subarctic 

Winter  Profile 

Visibi 1 ity 

40 

30 

20 

15 

10 

Wi  nd  = 

O 

0.963535 

0.963535 

0.963587 

O. 963549 

0.963600 

Wind  = 

5 

0.963548 

0.963644 

0. 963652 

0.963692 

0.963209 

Wind  = 

IO 

0.963594 

0.963567 

0.963533 

0.963565 

0.965684 

Wind  = 

15 

0.963558 

0. 963605 

O. 963507 

0.963501 

0. 956474 

Wind  = 

20 

0.963587 

0.963634 

0.963474 

0.963308 

0.956391 

U.S.  Standard  Atmosphere  Profile 

Vi sibi 1 i ty 

40 

30 

20 

15 

10 

Wind 

=  0 

0.920811 

0.920821 

0.920850 

0.920789 

0.920866 

Wind 

=  5 

0.920841 

0.920853 

0.920746 

0.920921 

0.922365 

Wind 

=  10 

0.920800 

0.920763 

0.920945 

0.920259 

0.919745 

Wind 

=  15 

O. 920802 

0.920818 

0.920675 

0. 920402 

0. 923306 

Wind 

=  20 

0.920836 

0.920866 

0.920934 

0.920379 

0.914907 
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TABLE  A. 8 


TABLE  OF  B— COEFFICIENTS  FOR  SIX  LOWTRAN  PROFILES  AT  3.5  um. 


Tropical  Profile 

Visibility  40  30  20  15  lO 


Wind  =  O  -0.078300  -0.079883  -0.083265  -0.086652  -0.093444 
Wind  =  5  -0.132564  -0.154864  -0.199849  -0.245554  -0.337061 
Wind  «  lO  -0.143678  -0.170376  -0.223794  -0.277764  -0.384270 
Wind  *  15  -0.146815  -0.174519  -0.230264  -0.286954  -0.401796 
Wind  =  20  -0.148416  -0.176768  -0.233486  -0.291292  -0.404436 


hid latitude  Summer  Profile 

Visibility  40  30  20  15  lO 


Wind  =  O  -0.064275  -0.065945  -0.069313  -0.072711  -0.079494 
Wind  =  5  -0.118579  -0.140899  -0.185840  -0.231112  -0.322254 
Wind  =  lO  -0.129858  -0.156395  -0.209959  -0.263693  -0.371488 
Wind  *  15  -0.132919  -0.160663  -0.216479  -0.273163  -0.384626 
Wind  -  20  -0.134394  -0.162714  -0.219772  -0.277832  -0.386820 


hidlatitude  Winter  Profile 

Visibility  40  30  20  15  lO 


Wind  =  O  -0.028635  -0.030311  -0.033691  -0.037106  -0.043922 
Wind  *  5  -0.082925  -0.105237  -0.150164  -0.195458  -0.286223 
Wind  =  10  -0.094293  -0.120931  -0.174555  -0.228508  -0.338264 
Wind  =  15  -0.097416  -0.125235  -0.181256  -0.238064  -0.352103 
Wind  =  20  -0.098874  -0.127235  -0.184428  -0.242007  -0.359376 


Page  1  of  2 


TABLE  A. 8 


TABLE  OF  B— COEFFICIENTS  FOR  SIX  LOWTRAN  PROFILES  AT  3.5  um. 


Subarctic  Summer  Profile 

Visibility  40  30  20  15  lO 


Wind  *  O  -0.048863  -0.050536  -0.053881  -0.057260  -0.064021 
Wind  ■  5  -0.103225  -0.125554  -0.170457  -0.215911  -0.305416 
Wind  *  lO  -0.114354  -0.140914  -0.194472  -0.248144  -0.359279 
Wind  =  15  -0.117388  -0.145095  -0.201005  -0.257037  -0.370669 
Wind  =  20  -0.118852  -0.147091  -0.203899  -0.261375  -0.372554 


Subarctic  Winter  Profile 

Visibility  40  30  20  15  lO 


Wind  *  O  -0.017768  -0.019474  -0.022912  -0.026375  -0.033305 
Wind  *  5  -0.071933  -0.094238  -0.139144  -0.184483  -0.274853 
Wind  =  lO  -0.083740  -0.110510  -0.164416  -0.218868  -0.328210 
Wind  *  15  — O. 087009  -0.115030  -0.1 71452  -0.228422  -0.340038 
Wind  =  20  -0.088548  -0.117162  -0.174745  -0.233146  -0.346934 


U.S.  Standard  Atmosphere  Profile 
Visibility  40  30  20  15  lO 


Wind  =  O  -0.036796  -0.038290  -0.041297  -0.044335  -0.050403 
Wind  *  5  -0.091927  -0.004360  -0.159548  -0.205238  -0.296B48 
Wind  =  10  -0.100272  -0.125884  -0.177567  -0.229439  -0.333287 
Wind  =  15  -0.102441  -0.128875  -0.182099  -0.235793  -0.344040 
Wind  =  20  -0.103449  -0.130276  -0.184322  -0.238687  -0.345855 


APPENDIX  B 


COMPUTER  CODING  TO  INSTITUTE 
ATMOSPHERIC  TRANSMITTANCE  CALCULATION 


The  fol lowing  FORTRAN  code  is  designed  to  institute  the 
calculation  of  atmospheric  transmittance  on  an  incremental 
range  movement  basis.  As  the  missile  closes  its  target (s) 
periodic  recalculation  o-f  acquisition  probabilities  based  on 
updated  range  and  atmospheric  transmittance  are  necessary  to 
realistically  model  an  engagement.  Listed  here  are  the 
subroutines  necessary  to  calculate  the  transmi ttance  based 
on  the  value  o-f  the  current  range  value  in  the  model. 
Although  these  subroutines  are  written  in  FORTRAN,  they  are 
simple  and  easily  convertable  to  other  languages.  Comments 
regarding  user  alternatives  are  contained  therein. 


C  First  thing  to  do  is  allow  the  user  to  select  which 
C  Atmospheric  Pro-file  is  desired.  Then  pass  this  to  the 
C  ATMOSP  subroutine  to  calculate  specific  coefficients. 

C  This  input  may  be  converted  to  disk-read  input  or 
C  NAMELIST  input  if  interactiveness  is  not  desired. 

INTEGER  ATPROF 

READ (111,*)  ATPROF 

111  FORMAT < '  Which  Atmospheric  Profile  is  desired  ?'// 

1 '  1.  Tropical  7 

2'  2.  Midlatitude  Summer'/ 

3'  3.  Midlatitude  Winter'/ 

4'  4.  Subarctic  Summer'/ 

5'  5.  Subarctic  Winter*/ 

6'  6.  U.S.  Standard  Atmosphere ', I 1 ) 


CALL  ATMOSP (ATPROF) 


C  The  values  of  Wind  and  Visibility  must  already  have 

C  been  input  to  the  program  prior  to  the  above  CALL. 

C  ATMOSP  calculates  the  ‘A'  and  *B'  coefficients  which 

C  are  necessary  to  refer  to  the  TRANSM  function. 


SUBROUT  I NE  ATMOSP  (ID 

DIMENSION  A1 (6) ,A2<6) ,A3(6> ,A4(6> ,C1 (6) ,C2(6> , 

1  D1 (6) ,D2(6) , SIGMA (6) 

COMMON  ACOEFF , BCOEFF , V I S I B , VW IND 

OPEN ( UN I T*9, FILE* 'ATMPROF.DAT' , STATUS* 'OLD ' ) 

DO  1  1*1,6 

READ (9, *)  A1(I) 

READ(9,»)  A2 ( I ) 

READ (9,*)  A3 ( I ) 

READ (9,*)  A4 ( I ) 

READ (9,*)  Cl (I) 

READ (9,*)  C2 ( I ) 

READ (9,*)  D1(I) 

READ (9,*)  D2(I) 

READ (9, *)  SIGMA (I ) 

1  CONTINUE 

CLOSE ( UN I T=9, STATUS*  KEEP ' ) 

ACOEFF  *  SIGMA (II) 

VCOEFF  *  Al(II)  +  A2 (ID *VISIB  +  A3 (II) *VISIB**2. 

#  +  A4 (II) *VISIB**3. 

WCOEFF  *  1.  ♦  D1(II)*(1.-  EXP (D2 (II) *VWIND) ) 

IF  (VWIND.LE.O. 1)  THEN 

VWCOEF=l . 0  !/*  No  correlation  if 

GO  TO  2  !/*  wind  near  zero. 

END  IF 

VWCQEF  *  Cl (II)  *  VISIB  **  C2 (I I ) 

2  BCOEFF  =  VCOEFF  *  WCOEFF  *  VWCOEF 

RETURN 

END 

C  The  two  necessary  coefficients  (A  or  o,  B  or  «)  are 
C  computed  and  global  for  referencing  by  TRANSM  for 
C  transmittance  calculation  based  on  current  range 
C  values.  Now  the  function  to  perform  equation  4.3. 
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FUNCTION  TRANSM  (RNG, AT) 

C  The  current  range  value  is  passed  in  as  RNG ,  and  AT  is 
C  the  Atmospheric  transmittance  value  passed  back  out. 

COMMON  RO , ACOEFF , BCOEFF 

AT  =  ACOEFF  *  EXP (BCOEFF *RNG) 

RETURN 

END 


The  file  ‘ATMPROF.DAT'  referred  to  in  subroutine  ATMOSP 
above  is  an  ASCII  data  file  which  contains  all  of  the  curve 
fitting  coefficients  listed  in  table  4.8.  It  is  easily 
constructed  with  an  editor,  keeping  in  mind  that  if  a  name 
change  is  necessary,  it  must  be  performed  in  both  places. 
The  format  is  shown  here  and  may  be  completed  by  referring 
to  table  4.8.  Ensure  that  no  comments,  blank  lines  or 
leading  spaces  occur  when  constructing  the  file. 


-O. 139154 

!/* 

Tropical 

A1 

6.2887E-03 

!/* 

A2 

-2.2354E-04 

!/* 

A3 

2.6091E-06 

!/* 

A4 

13.74105 

!/* 

Cl 

-0.716806 

!  /* 

C2 

0.895 

!  /* 

D1 

-0.30 

!/* 

D2 

0. 9383 

!/* 

9 

-0. 1249852 

!/* 

Mi d 1  at  Summer 

A1 

etc.  for  others  in  the  order  specified 
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